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Abstract 

The jet bundle description of time-dependent mechanics is revisited. The constraint al- 
gorithm for singular Lagrangians is discussed and an exhaustive description of the constraint 
functions is given. By means of auxiliary connections we give a basis of constraint functions in 
the Lagrangian and Hamiltonian sides. An additional description of constraints is also given 
considering at the same time compatibility, stability and second order condition problems. Fi- 
nally, a classification of the constraints in first and second class is obtained using a cosymplectic 
geometry setting. Using the second class constraints, a Dirac bracket is introduced, extending 
the well-known construction by Dirac. 
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1 Introduction 



For many years the problem of quantizing singular Lagrangians has been the object of much at- 
tention. Indeed, if we have a regular Lagrangian £:TQ — > K defined on the velocity space TQ 
of a configuration manifold Q, we get a nice Hamiltonian description on the phase space T*Q: the 
Lagrangian energy is transported from TQ to T*Q by the Legendre transformation and, using the 
canonical symplectic structure on T*Q we obtain the corresponding Hamilton equations. Next, we 
can use the standard quantization rules for the canonical Poisson bracket on the cotangent bundle 



When the Lagrangian is singular, however, (i.e. the Hessian matrix of £ with respect to the 
velocities is nonregular) we do not have a nice Hamiltonian description. In fact, not all the momenta 
are available, so we have some primary constraints defining a submanifold of T*Q. Moreover, given 
an initial data and its possible evolution, we have to ensure that it would be again admissible (the 
tangency condition) . Dirac JO] solved the problem by developing a constraint algorithm which gives 
(in the favourable cases) a final constraint submanifold where a solution of the dynamics exists. In 
addition, Dirac classified the constraints into two categories: first and second class constraints. The 
second class constraints were then used to modify the canonical Poisson bracket to obtain a new 
one (now called the Dirac bracket) which allows us to use the quantization rules. In this approach, 
an accurate description of the constraint functions plays a crucial role. Dirac's main aim was to 
apply this procedure to Field theories; indeed, many field theories (for instance, electromagnetism) 
are singular. 

After Dirac, a lot of work was done in order to geometrize his algorithm. The first important 



step was the work by Gotay et al [19], and its application to the Lagrangian formalism [20, 21]. 
Extension to field theories was always the main aim. Other algorithms were given later, in order 
to finding consistent solutions of the dynamical equations in the Lagrangian formalism of singular 
systems (including the second-order problem) (3|, and afterwards, new geometric algorithms 
were developed to be applied both in the Hamiltonian and the Lagrangian formalisms [BL |23|, |35|, 



,11 



However, the "Hamiltonian" description of field theories, termed the multisymplectic approach, 
is the natural extension of time-dependent mechanics. Therefore, if we wish to understand the 
constraint algorithm for field theories in a covariant formalism, the first step would be to extend 
the Dirac procedure to time-dependent Lagrangian systems. Some work was provided in 0, 10, [l7], 



24, 31, B3, 34, pq, Eq]. In this paper, we present a complete covariant description for this kind 
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of system, which generalizes the results of some of the above-mentioned references. Let us explain 
this formalism: 

We consider a configuration fibred manifold tt : E — * B, where B is an oriented 1-dimensional 
manifold (M or S 1 ) with volume form r/. The Lagrangian density is C = £rj, where £ is a function 
£ : J 1 E — ► R. Here J l E denotes the 1-jet prolongation of E which is called the evolution phase 
space of the system. If the Lagrangian is regular, then the dynamics is provided by the Euler- 
Lagrange vector field Xc, which is the Reeb vector field of the cosymplectic structure (FtciV) on 
J x i?, where flc is the Poincare-Cartan two- form and 77 is the pull-back to ^E of the volume form 
on B: 

i(X c )n c = 0,i(X c )ri = l. (1) 

If C is not regular, then {Vic, v) is 110 longer cosymplectic so that Eq. (|l|) has no solution in general 
(in addition, if such a solution exists it is not necessarily a second order differential equation). As- 
suming some weak regularity (constancy of the ranks of some distributions), a constraint algorithm 
can be developed in such a way that, in the better cases, we obtain a sequence of submanifolds 
which ends in some final constraint submanifold Mf of J^E. This is done in Section 3. 

The constraint functions defining these submanifolds were carefully determined in JlO| when 
E is trivially fibered, say tt = pr\ : E = R x Q — ► R, in which case J X E = R x TQ. Using a 
convenient decomposition of the Poincare-Cartan 2-form Vic based in the canonical connection in 
that trivial fibration, the authors have characterized the constraint functions. In this paper, we use 
a decomposition of flc based on an arbitrary auxiliar connection in J^E" — > B. The connection 
allows us to choose a basis of constraints, and could be understood as a choice of reference. 

We also examine the Hamiltonian setting (Section 4), and develop the corresponding constraint 
algorithm. We give a characterization of the constraints which can be connected to the Lagrangian 
ones by means of the Legendre transformation TL. This can be done in the case of almost-regular 
Lagrangians. The characterization obtained here is based again on the choice of a connection in the 
bundle V — ► B, where V = TL{J l E) is a submanifold (the primary constraint submanifold) of 
the restricted momentum dual bundle J l *E. As in the Lagrangian setting, the connection enables 
us to choose a basis of constraints. In addition, the Hamiltonian dynamics is obtained and related 
with the Lagrangian dynamics. We also consider the second order differential equation problem, 
since the Euler-Lagrange equations are of second order. 

Let us say that both algorithms can be considered in abstract as particular cases of a general 
algorithm developed in Section 2 for the case of a precosymplectic fiber manifold F — > B. 

In Section 5 we consider the problem of finding the dynamics, but imposing from the very 
beginning that this dynamics should simultaneously satisfy the second-order differential equation 
condition. We characterize the new constraints and relate them to the above approach. The results 



of this Section are a generalization to the time-dependent case of those given in |J] and [38| for 
autonomous systems. 

In Section 6, we classify the Hamiltonian constraints into first and second class, and define some 
geometrical projectors which allows us to introduce a Poisson bracket (the Dirac bracket) that has 
properties similar to the Dirac bracket in the autonomous case: second class constraints become 
Casimir functions, and the evolution of an observable is given by the bracket with a suitable 
Hamiltonian function. Note that now we need, in addition, a suitable vector field in order to 
describe the evolution (the bracket is not enough to do this task). 



Finally, we include two examples in Section 7 to illustrate how the procedure for chosing con- 
straints works. 
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The paper contains also three auxiliary Sections with the purpose of making it selfcontained. 

It is very important to point out that, throughout the paper, an auxiliar connection is used 
for constructing different geometrical structures. In fact, this technique (the use of a connection) 
was used for the first time in ||, in order to obtain (global) Hamiltonian functions, and afterwards 
applied both in the Lagrangian and Hamiltonian formalisms for this and other purposes (see |l3| , 

We trust that the results contained in the present paper will provide new insights into the 
problem of characterizing the constraints in classical field theories, and, consequently, in their 
quantization procedure. 

Manifolds are real, paracompact, connected and C°°. Maps are C°°. Sum over crossed repeated 
indices is understood. 



2 The general case 

2.1 Statement of the problem and solution 

Let k: F — > B be a fiber manifold, where dim F = 2N + 1 and the fibers have even dimension 2N 
(B can be M or S 1 ). Let w G Q 1 {B) be a volume form, and denote n = k*xu. 

Let V be a connection in k: F — > B; that is, V is a /-c-semibasic 1-form on F with values in TF 
such that V*/3 = (3, for every At-semibasic form j3 G i? 1 (i ? ). As is known, a connection V always 
exists, and it defines a horizontal subbundle H(V) C TF, such that TF = H(V) © V(k), where 
V(ft) is the k- vertical subbundle. If x G F, then H X (V) = Im V x . If y G X(F) is the vector field 
spanning the horizontal subbundle H(V) such that i(y)rj = 1. Then we can write V = r] ® 3^- 

For every X <E X(F), 

i(X)V = i(X)(rj (8 y) = (i(X)rj)y = X H e X(F) 

is an horizontal vector field; that is, a section of H(V) — * F. X H is the horizontal component of 
X, and we write X = X H + X v , where X v = X - X H £ X V(k) (F) (it is a k- vertical vector field). 

Furthermore, if a E f2 r (F), then 

i(V)a = i(rj <g> y)a = r) A i{y)a = a H . 

We have that a = a H + a v , where a v = a — a H is a k- vertical r-form with respect to the connection 
V, that is, it vanishes under the action of the horizontal vector fields associated with the connection 
V, and in particular i(y)a v = 0. Moreover, if X\, . . . ,X r G X(F) are k- vertical vector fields, then 
a H (Xx,...,X r ) =0. 

The problem we wish to solve can be posed in the following way: 

Statement 1 Given k:F—>B, and n as above, let Q, G Q 2 (F). We wish to find a submanifold 
C <— » F and a section X: C — > TqF, verifying that X x G T X C, for every x G C , and such that the 
following equations hold: 



[i(X)n] x = , [i{X)i 1 ] x = l , for every xGC 



(2) 
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A section X: C — > TcF is a "vector field on F with support on C" (we will denote the set of 
these vector fields by X(F, C)). 

This is the problem given by (F, f2, rj). Now we are going to state the above system of equations 
in an equivalent way. 

Theorem 1 Consider (F,$l,rj). Let V be a connection in k:F — > B, and consider the induced 
splitting ofVt. Then, the couple (C, X) (where C F is a submanifold, and X G 3C(F,C) is such 
that X x G T X C, for every x G C) is a solution for the problem stated by equations (^j if, and only 
if, the following system holds: 

[i(X)n V ] x = -[i(y)n] x , [i(X)rj] x = l ; for every xGC. (3) 



( Proof ) We have that 

n = n H + n v = v a [i(y)n] + n v 

and then, if X G X(F), 

i(x)n = d(x) v )(i(y)n) - r,[i{x) i{y)n] + i(x)n v . (4) 

If (C, X) is a solution of equations @, from (Q) we have that 

= [i(y)n] x - r, x [i{X) i(y)n] x + [i(X)n v ] x , for every x G C 
but [i(X) i(y)Q] x = 0, because [z(X)f2]a; = 0, and the equations (|3|) hold. 

(<£=) If (C, X) is a solution of equations (||) then, from @ we have that 

[i(X)n] x = -n x [i(X) i(y)n\ x = Vx [i(X) i(X)n v ] x = , for every x G C . 



From the equivalence of equations (||) (which are independent of the choice of the connection 
V) and (H), we have: 

Corollary 1 The couple (C,X), solution of equations (^), is independent of the connection V. 

The case we will focus on is: 

Assumption 1 Q, G Q 2 {F) and n G Q l {F) are closed forms, (andV = rj®y is a fixed connection 
in the fiber manifold k: F — > B). 

Then, we will denote 

j = i(y)n , uj = n v = n-n aj . 

(Notice that i(y)j = 0, and i(y)uj = 0). Finally, with this notation, equations (||) become 

[i(X)u] x = -7a, , [i(X)rj\ x = 1 ; for every x G C . (5) 
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Now, introducing the map (see Appendix ^]) 

(x,v) i-» i(v)uj x + (i(v)r] x )r) x 
we can state the following fundamental result: 

Theorem 2 The necessary and sufficient condition for the existence of a submanifold C > F 
and X G X(F,C), with X x G T X C, for every x G C, suc/i i/tai i/ie equations (||) (or equivalently 
equations ^)) hold, is that 

Vx ~ lx G t-(^,^)(TxC) ; for every x^C (6) 

or, iu/iai is equivalent 

{Vx ~ lx,^x~C) =0 ; /or every x G C 
w/iere b^^j = b^T.*- and T^C = (T^C) 1 - := [b^^C^C)] . 

( Proof) {=$■) It is an immediate consequence of Theorem |] and the definition of „>„.)• 

(<^=) Suppose that a submanifold C ^ F and a vector field X G X(F, C) tangent to C 
exist, such that @ holds. This means that, for every x G C, 

% - lx = \)( Vx ^ x) (X x ) = [i{X)u\ x + [i{X)n] x n x . (7) 

Contracting both members with y we obtain 

[i(y)vh - H(yh]x = [i(y)i(x)uj] x + [i(y)r,i(x) v } x 

and, as i(y)r/ = 1, ^(3^7 = 0, and i(y)u = 0, we conclude that [i(X)rj\ x = 1. Therefore, from (ffi) 
we have that [^(Xja;]^ = — j x , and the result follows from Theorem |]. ■ 

Remark 1 Every X G 3C(F,C) can be extended to define a vector field X G 3C(F). Therefore, 
although we are really interested in finding vector fields on F at support on C, from now on we 
will suppose that the vector fields are defined everywhere in F. 

2.2 The constraint algorithm 

Now we will apply the last result in order to solve the problem stated above; that is, to find a 
submanifold C *• F, transverse to the projection k:F—>B, and a vector field X G 3£(F), such 
that 

1. [i(X)n\ x = 0, and [i(X)rj\ x = 1, for every x G C. 

2. X\c is tangent to C. 

The procedure is algorithmic, and produces a sequence of subsets {Cj} of F. Then, we will assume 
that: 

Assumption 2 Every subset C{ of this sequence is a regular submanifold of F, and its natural 
injection is an embedding. 
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Thus, we consider the submanifold C\ ^ F where a solution exists; that is, 

d = {x G F | 3X G T X F : i{X x )Sl x = 0, i(X x )n x = 1} . 

Then there is a vector field X G X(F) such that = 0, [i(-X")f7]|ci = 1- But in general X 

is not tangent to C\. We can consider therefore the submanifold 

C 2 = {x G Ci I 3X G T^d : i(X x )Q x = 0, ^X,)^ = 1} . 

Then there will be a vector field X G 3t(F) tangent to C\ such that [i(X)0]|c 2 = 0, [i(^)^]|c 2 = 1- 
Again, X may not be tangent to C^- Following this process, we obtain a sequence of constrained 
submanifolds 

...'bc/CAcx^Co^F. (8) 
For every i > 1, Cj is called the zth constraint submanifold. 

This procedure will be called the constraint algorithm, and we have three possibilities: 

• There exists an integer k > such that = 0. This means that the equations are not 
consistent; that is, they have no solution. 

• There exists an integer k > such that ^ 0, but dimCfc = 0. In this case, there is no 
dynamics. consists of isolated points, and the solution of the equations is X = 0. 

• There exists an integer k > such that Ck+i = Ck = Cf, and dimC/ > 0. In such a case, 
there exists a vector field X G X(-F), tangent to Cf, such that 

[i(x)n]\ Cf =o, [i(x) v ]\ Cf = i 

In this case, the manifold Cf is called the fmai constraint submanifold. This is the situation 
which is of concern to us. 

Next we wish to give an intrinsic characterization of the constraints which define the constraint 
submanifolds d. In order to do this, following the same pattern as in Theorem g, we conclude 
that: 

Proposition 1 Every submanifold Ci (i > 1) in the sequence ^) can be defined as 

C i = {xeC i . 1 | (r^-T^T^a-i) =0} . 

Therefore, if the distribution T^Cj_i has constant rank, and {z'f , . . . , z\ % ^} is a set of vector 
fields spanning locally this distribution, then Ci, as a submanifold o/Cj_i, is defined locally as the 
zero set of the functions xf £ C°°(F) given by 

X f=i(zf-V)(r,-i) , (i = l,...,r). 
These functions are called ith-generation constraints. 
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2.3 Stability of vector fields solution. New characterization of constraints 

The solutions of the equations @ are now analyzed. Thus, consider the problem posed in State- 
ment |l[ Once we have found the submanifold C\ where a solution exists, we impose the tangency 
condition for the vector field solutions. This enables us to obtain a new characterization of the con- 
straints defining the submanifolds of the sequence (8|). To do this we need the following additional 
hypothesis (and the results stated in Appendices [A and |C]): 

Assumption 3 The distributions T^Co = T F = ker a; n ker rj, and T^ Co D TCi have constant 
rank. 



From the first part of Lemma 13 we deduce that these conditions are independent of the con- 
nection V. 

If X G X(F) verifies that 

[i(A>]| Cl = - 7 | Cl , [i(X) V )\ Cl = l (9) 

then, from now on, we can assume that i(X)r] = 1 everywhere in F. In fact, since if i{X)n ^ 1, we 
take X' = X + [l-i(X)rj\y, and it is clear that X'\ Cl = X\ Cl , &ndi(X')n = i(X)r] + l-i(X)n = 1. 
Therefore: 

Proposition 2 Let X be a vector field on F such that equations ^) hold, and verifying the above 
condition. 

1. If Z,Y G ker uj n ker rj, such that Z is tangent to C\, then 

i([x,z\ - [Uy)i(zyy]y)u\ Cl -v[i([x,z\ - [LOWWMk = &[i(zh]\ Cl 

and, as a consequence, 

{[X, Z) - [L(y) i(Zh)}y + Y} x G T^Ci , for every x G C x . 

2. If Z' G 3£(Cb) such that Z'\ x G T^Ci, for every x G C\, then there exist locally Z,Y G 
ker uj fl ker -q, such that Z is tangent to C\, and Z'\c 1 = {[X, Z] — [L{y) i(Z)j]y + ^}|ci • 



( Proof) 



1. It is clear that Y x G T X C±, since Y G ker oj n ker n, hence \>( VxyLUx ^(Xi)(Y x ) = 0, for X\ G 
T Zl Ci. 

For [A,Z] - [L(y)i(Z)-f)]y, we have 

i([A,z] - uy)%{z)i)\y)n = l(x) i(z) v - i{z) l(x)t7 - [l(30 »(z) 7 ] 

= -L(y)»(z) 7 1 J 

and 

Z] - i(^)7]y)w = ifl*, Z])w - [L(y) i(Z)i\) i(y)u = i([X, Z])uj . (11) 
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However, from the second equality, for every X' G X(Co), and bearing in mind that Z G 
ker u n ker r\ and that u = ft, — rj A 7, we have 

i(X') i([X, Z])u = i(X')[UX) i(Z) - i(Z) L(X)]uj = i(X')[- i(Z) i(X)du - i{Z)di{X)uj] 
= i(X')[-i(Z)i(X)dn-i(Z)i(X)(riAdj)-i(Z)di(X)u;} 
= - i(X') i(Z){[i{X)ri\d-y — rj A i{X)d-f + di{X)u} 
= -[i(X)n]i(X')i(Z)6rf- [i(X')ri\[i(Z)i(X)di\ -i(X')i(Z)di(X)u 
= -[t(X)n}i(X')i(Z)d 1 -[i(X')r ] }[i(Z)i(X)d 1 ) - 

i(X') L(Z) i(X)uj + i(X')d %{Z) i{X)u 
= -[t(X)n}i{X')i(Z)d 1 -[i(X')r ] }[i(Z)i(X)d 1 ) - 

L(Z) i(X') i(X)u + i([Z, X'}) i{X)uj . 

On Ci, using @, and taking into account that -h{Z)[i(X') i(X)u}\ Cl = [h{Z) i(X')j}\ Cl 
(since Z is tangent to Ci), we obtain 

i(X')i([X,Z])u\ Cl = -i{X') l {Z)d 1 \c 1 -[i{X')ri}{i{Z)i{X)d 1 }\ Cl - 

L(Z) i(X') %{X)u\ Cl + i([Z, X'}) i(X)u\ Cl 

= -i(x')L(z) 7 | Cl +i(x')di(z)~f\ Cl - [i(x') v ][i(z)i(x)di\\ Cl + 

UZ)i(X'h\ Cl - i ([Z,X'Mc 1 
= l {[Z,X']) 1 \ Cl +UX')i{Z) 1 \ Cl - 

[i(X')r ? ][ i (Z)f(X)d 7 ]|c 1 -i([^^ , ])7lc 1 
= UX^HZhlc.-iiiX'MiWUX^lc, ■ (12) 

Furthermore, using that 7 = i(y)Q and that £1 is closed, 

i(z)i(x)d-f\ Cl = i(z)i(x)L(.y)n\ Cl 

= L(y)i(z) i (x)n\c 1 +z([y,x})i(z)n\ Cl -i([y,z}) i (x)n\ Cl 

= L(y)i(Z)j\ Cl , (13) 

since i(X)£l\d = (because X is a solution), i(Z)£l\c± = (by Lemma |l8|, Appendix |c|), 
and i(Z) i(X)U = i(Z) i{X)uJ + j(Z) 7 = i(Z) 7 . Therefore, from flu]), @ and (|J), 

Moreover, from (|l0|), 

i(X'){^([X,Z] - [L(^)i(^)7]^)r?]}|c a = -[i(X')r/]L(^)i(Z) 7 )|c 1 • 
Hence we conclude 

i(x') [i([x, z] - [L(y) i(z) 7 ]y)u - v [i{[x, z\ - [L(y) i(z) 7 ]y) v ]] \ Cl = 

= [L(X')i(ZM Cl = [i(X')d[i(Z) 1 ]] \ Cl , 

for every X' G X(Co), so 

i([X,Z] - Uy)i(Zh]yU Cl -v[i([X,Z] - [L(y) i(Z)j}y)r]]\c\ = d[i(Z) 7 ]| Cl . 

Finally, notice that T^Ci = {v x G T X C Q \ (b( Vxj u> x )( u x),v x ) = , Vu x G T^Ci}. That is, 
(i(v x )oj x - [i(v x )r] x ]r] x ,u x ) = 0, for every u x G T X C\. Taking v x = {[X, Z] - [L(y) i(Z)j]y + 
Y} x , we obtain that {%{v x )u) x - [i(v x )r] x ]rj x ,u x ) = [L(W) i(Z)j]\ x = 0, if U is a vector field 
tangent to Ci such that f/(x) = u x (notice that i(Z)(rj — j) = — i{Z)^j is a constraint function 
defining C\ in Co). Therefore, v x G T^Ci. 
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2. Consider a local basis {zf\ . . . , Z? x } of the distribution ker tinker 77, in an open set U C F. 
Then the 1-forms di(Z- )7|c? 1 , i = 1, . . . , ri generate the annihilator of T^nCiCi- Therefore 



{i{Z>)u-[i{Z')ri)]n}\ Cl =n&i{Z^ ) 7 ]| 



Ci 



where /* € C°°(C/ D Ci) are arbitrary functions, which can be extended to functions /* G 
C°°(*7). Thus, we consider Z = fzf\ then i(Z) 7 = /S(zf } ) 7 , and 

[di(z) 7 ]k* = f\ Cl [di(zl 0) Md + (dh\ Cl i(z(%\ Cl = f | Cl [di(zf ) 7 ]| Cl . 

Hence, comparing these last equations, we conclude that 

{i{Z')u - [i(Z') V )] V }\ Cl = [di(Z) 7 ]\ Cl (14) 
and in this way, for every z| , 

= {i{zP)(i(Z')u - [i{Z')ri)]ri)}\ Cl = L(zf ) i(Z) 7 \ Cl . (15) 

Now, 

l (Z)Uzf ) )l\c 1 = {i(Z)i(zf ) )& 1 + l (Z)&i(Z%}\c 1 = {i(Z)i(Z^ )d 7 + L(Z)*-(Z°>) 7 }h • 
Also, 

i([zf\z]) 7 | Cl = {UZ? ) ) l (Z) 1 -i{Z)UZ ( ? ) h}\c 1 ■ 

Therefore, 

i(Z)^(zf ) )d 7 | Cl = {^) L (zf ) 7 - L (Z)*(zf ) 7 }| Cl = 

= { L (zf ) i(Z) 7 - L(Z) t (zf ) ) 1 - i([zf\z]) 7 }\ Cl . 
Taking into account the first item of Lemma [l^ (Appendix |C|) , we have 

= i(Z)i(zf>)d 7 |c 1 =HzP)i(Z) 7 \ Cl -nz)i(zi%\ Cl - t ([zf\z)) 1 \ Cl . 

Hence, L(Z) i{zf >) 7 | Cl = L (^ 0) ) i(Z) 7 | Cl - K^f, ^bb • But L(^ 0) ) i(Z) 7 | Cl = 
(see (|l^)) and, using the second item of the Lemma [l9| and the characterization of C\, 
i([Z-°\ Z])7|c 1 = 0. Therefore, Z is tangent to C\. Finally, from ([[]), (|H]) and the first part 
of this Proposition, we conclude that 

{i{z* - [x, z] + [L(y) i(z) 7 ]y)cj - H(z' - [x, z] + [L(y) i(zh]y)n] v }\ Cl = 
= {i(z')u- Uz'^Mic, - {i([x, z] - [L(y) i(z)-y]y)u-[i{[x, z] - [L(y) i(z) 7 ]y) v )] v }\ Cl = 

= di(Z)-f\ Cl - di(Z)j\ Cl = 

and this means that 

Z' - [X, Z] + [L(30 i(Z)i\y G T^C = ker oj x n ker n x , Vx G Ci , 
i.e., = {[X, Z] - [L(y) i(^) 7 ]y + Y}\ Cl , with Y" G ker u D ker 77. ■ 



Using this Proposition we can construct a local system of generators for T C\. In fact: 



Corollary 2 
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Let X G X(F) be such that equations ^ hold, tq the rank of the distribution T^Co fl TC\, and 
{z[ ^ , . . . , Z$ , Z^ +1 , . . . , Zf® } a local basis for ker tinker rj (in a neighbourhood of a point x G C\), 
where {Z± }i=i ... ro are tangent to C\. Then T^Ci is locally generated by the vector fields 

{zf\ z£U£k, ■ ■ • . 4% zf] - [ L (30 i{z\%]y} i=1 _ ro . 



Remark 2 From now on, we will denote by 3C (Cj) the set of vector fields of X(F, Cj) which span 
locally the distribution T^Ci. This set can be characterized as 

X ± (C i ) = {ZeX(F) | j*[i(Z)co-(i(Z) V )n]=0} 

where C% <— > F denotes the corresponding embedding. 



Since C2 = {x G C\ \ (r] x — ^ X ,T X C\) = 0} and i(Zj )(r] — 7)|c*i = , if we introduce the 
functions 



xf ■■= i ([x,z^} - Uy)i(z^)j}y)( v -j) 

we conclude that a basis of constraints for C2 is made by those functions which do not vanish on C\ . 
Bearing in mind that i([X, Z^])n\c 1 = 0, i(y)^f\c 1 = and i(y)rj = 1, they can be characterized 
as 

x p = - i ([x,zP}) 1 - (y[i(zi%]). (16) 

Next we will analyze the solutions of the dynamical equations. First we describe the set of 
solutions on C%. 

Proposition 3 The general solution of the equations ^) on the submanifold C\ <— > F where they 
are compatible is X^°> + Y^°> , where X^ G 3C(F) is a particular solution and Y^°> G X(F) is an 
arbitrary element o/ker Q n ker 77. 

( Proof) It is a consequence of the linearity of i(X)Q and i(X)rj on the vector field X. ■ 
Taking into account the above results, we can write these solutions in the form 

x = + y<°> = jr<°) + fzf + fzf (i = i, ... ,r ; k = r + ,n) . 

Now we must identify the points of C\ where there are solutions tangent to C\. Consider the set 
of independent constraints {x^} defining locally C\ as a closed submanifold of Co = F, which by 
Proposition [l] are described in the following way 

{i( z i )(V ~ 7)}i=l,...,r , {i{zf ] ){n - -i)} k=ro+ ^.^ ri . 
Then the tangency condition imposes that 

= L(X(°) + fzf + f k zf ) ) t {zf ) ){n- 1 )\ Cl 

- -L(x(°))i(zf ) 7 | Cl -f L (zf )^) 7 | Cl -f k L(Z^ ) i (zf) 1 \ Cl (17) 

= L(X(°) + rzf + / fe zf)^)(r/- 7 )| Cl 

= -L(^X(zW) 7 k -/ i L(zf)^) 7 | (7l -/ fc L(zf)^) 7 | Cl (18) 

where we have taken into account that, as a consequence of Lemma L(Z (0) )r7| Cl = 0, for every 
G X" L (Co). We will need the following result: 
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Lemma 1 The matrix (A^i) = (h^Z^ )i(zj^ )"/) (k,k f = ro + l,...,ri) is regular, for every 
x G C\ (in the corresponding neighbourhood). 

( Proof ) Let us suppose that the rows of this matrix are dependent for a point x G C\ . Then, 
u*zf{x x )(i{zf})i) = K , for all k', with / gl, 



Then, by Lemma 19, if Z- , i = 1, . . . , ro, are tangent to C\, we have that 

L(4 0) )^ 0) hk = *([zf ,zf)]) 7 | Cl +^ 0) )L(4 0) )7b 
= i(zf 0) )i(zi 0) )d7| Cl + *(zf )d(*(zf ) 7 )| Cl = L(^ 0) )i(zf ) 7 | Cl = 0, (19) 

and hence 

/i *Z J [ 0) (x)(i(^ 0) )(7)) = 0, for every i = 1, . . . , n . 

Therefore n\zf\ G T x Ci and hence n k {zf\ G T^C n T^Ci, for every x G Ci, but as 

{(zj ) x , . . . , (Zro ) x } is a basis of T^ Co n T x Ci, we must conclude that jjL k = 0, for every k = 
r + l, ■ 

Since L(4 0) )i(zf ) ) 7 | Cl = h{zf ] ) i(zf } ) 7 | Cl = 0, system © becomes 



= L(X(°)) *(zf ) 7 | Cl + /* L(4 0) ) *(zf) 7 | Cl 

which is a linear system on the coefficients f k . As a consequence of the last Lemma, this system 
has a unique solution that gives us the coefficients f k \ci- Then the general solution on C\ can be 
written as 

x = xw> + fzf + fzf = jstW + y« (i = 1, . . . , r ; fc = r + 1, . . . , n) 

where = X^ -* + f k Z^ and y^ 1 ) denote the determined and the undetermined parts of the 
solution, respectively. 

Furthermore, since 

L(^ 0) ) *(zf ) 7 | Cl = Wf \ Zf Mci + l{ zf ) L (zf ) 7 | Cl 

= -i(30<([zW,zW ])fi| Ci +i (zW )i{zf } )d 7 | Cl - *(zf )d[*(;y)*(zf )fi]| Cl = 
(where we have used that 7 = i(y)£l and Lemmas [j^ and |19|), system ( |T7| ) reduces to 

= L(Jf(°))i(4 0) )7l 



Ci 



and, if these equalities do not hold, they give the constraint which defines locally the submanifold 
Ci as a submanifold of C\. In fact, we have that 

L(x(°))[i(zf ) 7 ]| Cl = *(x(°))d[i(zf ) 7 ]| Cl 

= i(x(°)){,([x, zf ] - [loo *(4 0) )7]y)« 

- r,[i([x,z^]-[L(y)i(z^ h}y) v }}\ Cl 

= i(x(°))i([x,zf }p + Uy)i{zf mix^Hc, 

= -i([x,zp ])*(x(°)v| Cl +Uy)i{zf hid 

= i([X,z( 0) ]) 7 | Cl +L(^)i(zf ) 7 | Cl 
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where we employ the first item of Proposition || and i(X^)n\c = 1, i(X)r] = 1 (since X is a 
particular solution). Observe that this result is just the characterization of the 2?i(f-generation 
constraints given in (|l6|). 

This study has been done for the submanifold C%, but it can be extended to every level of 
the constraint algorithm. As in the case of the submanifold C\, we will assume that for every 
submanifold Cj (i > 0) of the sequence @, the distributions T^Cj and Tq Ci n TCj + i have 



constant rank (note that from the second part of Lemma 13, we deduce that the first of these 



conditions is independent of the connection V). Under the above conditions, we have 

Theorem 3 Consider the sequence {Ci} of submanifolds (^j. For every i > 1, let r(j_x),r.j be the 
ranks of the distributions T Q-i and T^C^i D TCj respectively, and consider a local basis for 
X" L (Cj_i) (in a neighbourhood of a point x G Ci) 

r z (i-l) ~(t-l) z (i-l) 7^- l H 

1^1 '•••' Z 'r (l _ 1) > z 'r (l _ 1) +l'---'^T I / 

where {Z^ 1 ^}j=i,...,r (i _ 1 ) c,re tangent to C{. Finally, let X = .X'f* -1 ) +y(' —1 ) 6e i/te general solution 
on the submanifold Ci, where X^ 1 ' and y(' t_1 ) denote the determined and the undetermined parts 
of the solutions, respectively. 

1. Every submanifold Ci + \ in this sequence can be defined (in Ci) as the zero set of the so-called 
(i + l)th-generation constraints, {x '} C C°°(F), which are characterized in the following 
equivalent ways: 

(a) For every 

x (l+1) =i(z^)(v-i). 

(b) For every j G {1, . . . ,r (i _ 1} }, 

x (m) =L(x ( i -i) )[ . (z (-D )(r? _ 7)] . 

fcj For every j G {1, . . . , r(j_i)}, 

x ^= i ([x^\zf- 1) ]) 1 + uyMzf- 1) h] . 

2. For every k G {r7j_i) + 1, . . . ,rj} the stability condition 
determines (partially or totally) y( i_1 ) (on Ci). 



3 Lagrangian formalism for mechanical systems 



In this section we wish to apply the above results in order to study the evolution equations of 
non-autonomous singular Lagrangian systems. By using the constraint algorithm we will obtain a 
final constraint submanifold defined by the so-called dynamical Lagrangian constraints. However, 
this problem shows special features since, in addition to solving the dynamical equations, we must 
find solutions satisfying the so-called Second-order condition. This condition will be studied in 
Sections |4~5| and 0. 
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3.1 General description of Lagrangian systems 

(For more details on the jet bundle description of non-autonomous Lagrangian systems see, for 
instance, Q, g3|, 0, 0, @, || and @). 

Let tt: E — > 1? be the configuration fibered manifold of a non-autonomous mechamical system, 
with volume form ot € J? 1 (2?), and E" is an (n + l)-dimensional manifold, tt . J^E — > E is the 
jet bundle of local sections of tt, which is called the evolution phase space of the system. The 
map 7T 1 = tt o tt 1 : J l E — > B defines another structure of differentiable fibered manifold. Finally, 
(t,q p ,v p ) (with p = 1, . . . , n) will denote natural local systems of coordinates in J 1 ^ adapted to 
the projection tt: E — > B, and such that the form r] := tt^vj can be locally written as r\ = dt. 

Note that J l E is an embedded submanifold of the tangent bundle TE of E. In fact, the map 
l : J X E -> TE defined by 

L{jt <t>) = <t>*to(V m (to)) 

for to £ B and a local section of tt : E — > B, is an embedding. Here, denotes the vector field 
on B characterized by the condition w(V m ) = 1. If we consider fibered coordinates (z,q p ,v p ) on 
J^E and (t,q p ,i,v p ) on TE then the local expression of i is 

i(t,q p ,v p ) = (t,q p ,l,v p )- 



The dynamical information is given by introducing a Lagrangian density which is a tt 1 -semibasic 
1-form on J^E. Then there is a function £ G C oc (J 1 E) such that £ = £rj, which is called the 
Lagrangian function associated with C and r\. The Poincare-Cartan 1 and 2-forms associated with 
the Lagrangian density C are defined using the vertical endomorphism V of the bundle J^E: 



e c :=i(V)C + Ce Q 1 {J l E) 
Then a Lagrangian system is a couple (J 1 ^, &c). 
In a natural chart in J l E we have that 

V = (dgP 



ft/ 



-d0 £ <E Q 2 (J 1 E) . 







and therefore 

n c 

Notice also that 
nAU n c = 



d£_A n p 
dvP a 1 



' d£ , 
k8vP 1 



£) dt , 



. -di/Adi/ 
ov p ov v 

( d 2 £ „ 

v 

\dq p dv u 



<} ~ 1 dq p A dq u + , v"dv p A dt+ 



&£ \ 



dv p dv h 



d£ 

dqP + dtdv p 



(20) 



J dq p A 



dt 



n(n + l) 

(-1) 2 n!det 



/ d 2 £ 
\dv p dv v 



dt A dg 1 A ... A dq n A dw 1 A ... A dt> r 



p,i/=i, 



As usual, we say that a Lagrangian function £ (and hence the corresponding Lagrangian system) 
is regular if its associated form ft£ has maximal rank, or what is equivalent, the couple (Qc,n) 

( d 2 £ \ 

is a cosymplectic structure. This is also equivalent to demanding that det I — — - — ) is different 

\ov p ov v J 

from zero at every point. There exists a more general notion of regularity in the recent geometrical 
approach based on Lepagean forms [28, 29, Q (see also f27|). For Lagrangian systems of order 1, 
both definitions coincide. 
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A variational problem can be posed from the Lagrangian density C, which is called the Hamilton 
principle of the Lagrangian formalism: the dynamical trajectories are canonical liftings of the 
sections of 7r which are critical for the functional L:F C (B, E) — > R defined by L(c/>) := f B (j 1 (j))*£, 
for every <p G T C (B, E) (T C (B, E) denotes the set of compact supported sections of n, and j 1 ^ is the 
canonical lifting of <f> to J l E). Therefore, it can be proved that the critical sections of the Hamilton 
variational principle are the integral curves of an holonomic vector field Xc G X(J 1 -E'), (also called 
a Second Order Differential Equation or SODE), satisfying that: 

i(x c )n c = o , i(x c )v^o 

The second equation means that the vector field Xc is ^-transverse. It is usual to write this 
condition in the form 

i(X c )v = 1 

and, locally, this is equivalent to fixing the parametrization of the integral curves, taking as param- 
eter the coordinate t. From now on we will follow this convention. 



Remark 3 Let <fi G T C (M, E) be a section such that is an integral curve of a vector field Xc 
which is a solution to the above equations. Then <p verifies the Euler-Lagrange equations: 



d£ 
~d(f 



d d£ 



dt dvf 







(for p = 1, 



and conversely. 



We can characterize SODE vector fields in J E in a more suitable way, by using the vertical 
endomorphism. In fact, as V G Q l {J 1 E)®T(J l E, V(7r 1 ))(g)r( J 1 ^, tt u TB), we can contract the last 
factor with elements o{T(J 1 E,ir u A 1 T*B), thus obtaining an element of Q 1 ( J l E) ®T( J l E, V(vr 1 )), 
in a natural way. Therefore, we define the canonical endomorphism J G Q 1 {J 1 E) ® T( J 1 ^, V(-7r 1 )) 
by making 

J = i{V)n 

whose local expression in a natural chart of J^E is 

J = i( V )r, = (dq p - v p dt) 8) A . 

Now, it can be proved that a vector field D G X(J 1 E) is a SODE if, and only if, the following 
conditions hold 

J(D) = and i{D)r, = 1 . 
So, in Lagrangian mechanics, we search for vector fields Xc G X(J 1 £') such that: 

1. They are solutions of the so-called dynamical Lagrangian equations (Q) 

i(X c )n c = , i(Xc)n = l. 

2. X c are holonomic (SODE); that is, J{X C ) = 0. 

As is well known, if {J l E, VLc) is a regular Lagrangian system, then there exists a unique Euler- 
Lagrange vector field Xc G X(J x i?) for this system (in fact, Xc is the Reeb's vector field of the 
cosymplectic structure (Oc, rj)). This result does not hold for non-regular Lagrangian systems, and 
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if (J 1 -^, £lc) is singular, (Qc,v) is no longer cosymplectic, and the equations (|T]) have no solution 
in general. Even if it exists, it will be neither unique nor a SODE. In the best of cases, these 
vector fields can exist only in a subset of points of J x i?, and the most interesting situation is when 
this subset is a sub manifold of J^E. Thus, the problem we want to solve is the following: to 
look for a submanifold S J X E, transverse to the projection n . J E — > B, and a vector field 
X c e X^EtS) such that 

1. [i(X c )ttc]y = 0, and [i(X c )rf\y = 1, for every y G S. 

2. Xc is a SODE on the points of S; that is, J{X£)\s = 0. 

3. Xc is tangent to S. 

Conditions 1 and 3 are called compatibility and stability or consistency conditions for the dynamical 
equations, respectively. The problem can then be solved by first considering only these conditions 
and afterwards adding the second one, or taking all of them simultaneously. 



3.2 Solving the dynamical Lagrangian equations. Lagrangian constraint algo- 
rithm and dynamical Lagrangian constraints 

First we consider the problem of finding a submanifold M <^-> J^E, transverse to the projection 
f 1 : J l E -> B, and a vector field X £ 3L{J l E,M) such that 

1. [i(X)n c ]y = 0, and [i(X)n]y = 1, for every y E M. 

2. X is tangent to M. 



To solve this problem, we apply the algorithmic procedure developed in Section 2.2, with the 



hypothesis there assumed. Following this process, we obtain a sequence of constrained submanifolds 



Mi 



Mi 



M n 



J L E. 



(21) 



For every i, Mi is called the ith Lagrangian dynamical constraint submanifold, and this procedure 
will be called the Lagrangian dynamical constraint algorithm. 

The only case that concerns us is when the algorithm ends by giving a submanifold Mf (with 
dim Mf > 0), which is called the final dynamical constraint submanifold for the Lagrangian prob- 
lem. In such a case, there exists a vector field Xc € X.(J l E), tangent to M/, such that 



[i(Xc)n c ]\ Mf =0, [i(Xc)r,]\ Mf 



(22) 



Next we wish to give an intrinsic characterization of the constraints which define the constraint 
submanifolds Afj. In order to do this, we apply the results stated in Section ^. First we take a 
connection V in the fibered manifold tt 1 : J l E — > B; that is, an horizontal subbundle of T(J l E) 
or, what is equivalent in this case, a ^-transverse vector field y G X(J 1 S), which can be selected 
such that i(y)r) = 1 holds. 



Remark 4 It is interesting to point out that this construction can be made starting from a con- 
nection V in tt: E — > B or, what is equivalent in this case, a 7r-transverse vector field 3^£? £ 3C(E). 
Then this connection induces another one in the fibered manifold tt 1 : J x E — > B, which is associated 
with the ^-transverse vector field 3^ = j £ X(J 1 £') (the canonical lifting of 3^; to J 1 ^). 
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Now, let 7£ G Q 1 (J l E) and uoc G l? 2 (J 1 i?) be the forms denned as 

lc ■= i{y)^c , u c := &c - V A ic ■ 

Thus, for every submanifold Mi and every y G Mi, TyMi is a vector subspace of Ty^E, and 
we can consider the orthogonal complement T^Mj with respect to the couple (rjy, (uc)y); that is 

Ty-M t = [b (r)s , (Wjc)s) (T s M i )]° = {y G TjJ 1 ^ | - Vy(i(Y)Vy)]h,M t = 0} . 

Note that T^Mq = ker n ker rjy. Denote by X (Mi) the set of sections of the vector bundle 
T^Mi -> Mj. Therefore, as a direct application of Proposition || and Theorem [3| to the present 
case, we have: 

Theorem 4 Let (J l E,Q.c) be a singular Lagrangian system, and consider the sequence {Mi} of 



submanifolds (21). For every % > 1, let r^_i),rj be the ranks of the distributions T Mj_i and 

T^.Mi-i n Ti 
pomi x G Mj j 



Tj^.Mi-i n TMj respectively, and consider a local basis for X (M{-%) («n a neighbourhood of a 



. (i-l) 57(1-1) jr(i-l)\ 

1^1 '•••>^'r fi _ 11 >^' m _n+l'---'^ri J 



• / r( i _ 1) >~r (i _ 1 )+l> 



where {Z- sj=x,...,r^_ 1 ^ are tangent to Mi. Finally, let X^ = X [ c + be the general solu- 

tion on the submanifold Mi, where 
parts of the solutions, respectively. 



tion on the submanifold Mi, where and denote the determined and the undetermined 



1. The submanifold M\ where the dynamical Lagrangian equations are compatible can be defined 
as the zero set of the so-called 1 si-generation dynamical Lagrangian constraints, {C^} C 
C°°(J 1 E), which are characterized as 

c (i) = z -(2(°))(ry - 7£ ) , for every G X^^E) . 

2. Every submanifold Mj+i in this sequence can be defined (in Mi) as the zero set of the so- 
called {i + l)th- generation dynamical Lagrangian constraints, {C (i+1) } C C°°(J 1 S), which are 
characterized in the following equivalent ways: 

(a) For every Z® G X ± (M i ), 

C (i+1) = - ( ^) )(r? _ 7£) . 

(b) For every j G {1, . . . ,r(j_!)}, 

C (^) = L( x2- 1) )[,(4 i - 1) )( ?? - 7£ )]. 

(c) For every j G {I, . . . ,r {i ^}, 

3. For every k G + 1, • • • , the stability condition 



C it+1) = l ([xt 1 \4~ 1) ]hc + Hyuzf-% c ] . 



o = l(4 <_1) + YtHazt^iv - ic)]\ Mi 



determines (partially or totally) X ' (on Mi). 
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4 Hamiltonian formalism for mechanical systems 

Now we consider the Hamiltonian formalism of a non-autonomous mechanical system associated 
with a singular Lagrangian and, in particular, the problem of finding solutions of the Hamilton 
equations. 

(For more details on the jet bundle description of the Hamiltonian formalism of non-autonomous 



mechanical systems see, for instance, p2fl, 113], [17], [18], [p2], [39] and El 



4.1 The momentum dual bundle and the Legendre map. Hyper-regular and 
almost regular systems 

In order to associate a Hamiltonian system with a non-autonomous Lagrangian system, we first 
need to introduce two new elements: a suitable dual bundle of J l E, and a Legendre map. 

In the jet bundle description of non- autonomous dynamical systems, there are several choices 
for the momentum bundle where the covariant Hamiltonian formalism takes place (see |lf| for 
a general review). In this work, following Q, first we take the bundle T*E, which is called the 
extended momentum dual bundle associated with n:E^B. The natural projections are denoted 
by <7 . T*E — ► E and a 1 : T*E — ► B. Then, if AqT*E denotes the bundle of 7r-semibasic 1-forms in 
E, we define the bundle 

J U E := T*E/A T*E ~ V*(vr) 

which is called the restricted momentum dual bundle associated with tt:E — > B. We denote the 
natural projections by r 1 : J E — > E, and r 1 := tt o t :J E — > 5. In addition, we have the 
canonical projection //: T*i? — > J X *E. 

Natural coordinates in T*E and J 1 *^ (adapted to the bundle structures) will be denoted by 
(t,q p ,p,p p ) and (t,q p ,p p ), respectively. Then we also write fj := f l *w = dt. 

The construction of the Legendre map can be made following different but equivalent ways. 
For instance, we can follow a procedure similar to the case of autonomous mechanics, and define 
this map as a "fiber derivative" of the Lagrangian density (see, for instance, Jll|). Alternatively, 
we can define it using the theory of affine dual bundles (as in || and |43j). Another more simple 
way consists in constructing this map by means of the Poincare-Cartan form as follows: taking into 
account that 0£ can be seen as a 1-form on ^E along the projection n 1 : J E — ► E, the extended 
Legendre map associated with a Lagrangian density C is the C°°-map T L: J l E — > T*E defined by 

[TC(y)](Z) = {Qc)y(Z) (23) 

where y € J 1 ^, Z € T^i^E, and Z £ TyJ l E is such that irly(Z u ) = Z u . Therefore, the restricted 
Legendre map associated with a Lagrangian density C is the C°°-map T C: J l E — > J^E defined by 

TL := p o TC 



In natural coordinates, the local expressions of these Legendre maps are 

„,p d£ 

(24) 



TLt = t , TLqP = qP , TC Pp = §§ , F£*p=£-vP§§ ! 



TL*t = t , FC*qP = qP , TC*p p = §§ 
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The cotangent bundle T*E is endowed with canonical forms: the Liouville forms O £ f2 1 (T*E) 
and f2 := — d© € f2 2 (T*E) (which is a symplectic form), whose local expressions are 

6 = pdt + p v dq v , n = -dp Adt- dp u A dq u . 

Using ( p0|) and (||), we obtain that 

fc*q = g c , TC*n = n c ■ (25) 



The matrix of the tangent map J-C* in a natural coordinate system is 



Id 

Id 

a 2 £ d 2 £ a 2 . 

dtdvP dq»dvP dVdvP 



Id (26) 

d 2 £ d 2 £ d 2 £ 



( <9 2 £ \ 

where the sub-matrix ( — — - — ) is the partial Hessian matrix of £. Thus, one deduces that (see 

\dv v dvP J 

ker^C* = kerFC* = ker £ n X V{ ^\j 1 E) (27) 

Note that the Lagrangian system {J E,Q.c) is regular, that is, (Qc,rj) is a cosymplectic structure 

(or equivalently the partial Hessian matrix I — — - — ) is regular everywhere in J 1 E) if and only if 

\ov v ovP J 

J-C is a local diffeomorphism (see also || for a different definition of this concept). Then, following 
a well-known terminology in mechanics we define: 



Definition 1 Let (J 1 ^,^) be a Lagrangian system. 



1. {J l E,Q.c) is said to be a regular or non-degenerate Lagrangian system if J-C is a local dif- 
feomorphism. 

As a particular case, (J 1 -^, ilc) is said to be a hyper-regular Lagrangian system if TC is a 
global diffeomorphism. 

2. Elsewhere {J l E,Vtc) is said to be a singular or degenerate Lagrangian system. 



Proposition 4 (See [15, Let {J l E,Q<c) be a hyper-regular Lagrangian system. Then: 

1. J-C{J l E) is a 1-codimensional embedded submanifold ofT*E which is transverse to the pro- 
jection \x. 

2. The map J-C is a diffeomorphism on its image, and the map \x restricted to J-C(J 1 E) is also 
a diffeomorphism. 



For dealing with singular Lagrangians we must assume minimal "regularity" conditions. Hence 
we introduce the following terminology: 

Definition 2 J7Z| , [?!]/ A singular Lagrangian system (J 1 ^, flc) is said to be almost regular if: 
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1. V := J-C(J E) is a regular submanifold ofT*E. 

(We will denote by jq:V T*E the corresponding embedding). 

2. J-C is a submersion onto its image (with connected fibers) . 

3. For every y G J^E, the fibers T X -1 \T C(y)) are connected submanifolds of J l *E. 

(This definition is slightly different from that in references [17|, and M])- 
Let {J l E,Q,c) be an almost regular Lagrangian system. Denote 

V := FC{J l E) . 

Then let Jq'.V > J^E be the canonical inclusion, and ^q:V — > V the restriction of the map [i. 
Finally, define the restriction mappings 

TLq: J x E -» V , F£ : J l E -» V . 

With these definitions, it follows that 

Proposition 5 (See Jl5], [3^]). Let (J l E,Qc) be an almost regular Lagrangian system. Then: 

1. For every y £ J X E, 

F£ ~\f£ (y)) = FCq\FCo{v)) . 

2. There exists a unique differentiable structure on V such that fio : V — > V is a diffeomorphism. 

3. V is a submanifold of J^*E, and jq:V J^E is an embedding. 

4- The restriction mapping J-£q : J^E — » V is a submersion with connected fibers. 



4.2 Hamiltonian system associated with an almost regular Lagrangian system 

If (J l E, Qjc) is an almost regular Lagrangian system, the submanifold jo'.V > J^E, is a fiber 
manifold over E (and B). The corresponding projections will be denoted Tq : V — > E and Tq : V — > B, 
satisfying that r 1 o j = Tq and f 1 o j = fg. We will denote rp := Tq*w = dt. 

Now, given the diffeomorphism ho = Hq , we define the Hamilton-Cartan forms 

@h = (Jo ° h )*9 ; O ho = -d9 ho = (Jo o h )*Q 

and, taking into account the commutativity of the diagram 

V ^ _ ji* E 

h H 
ji E ^ -^o . p 2° . 

from (^) we obtain that the following relations hold 

T£* e ho = G c , F£* n ho = n c . (28) 
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Then (J 1 * E,V,flh ) is the Hamiltonian system associated with the almost regular Lagrangian 
system {J l E, 

We can state a variational problem for ( J^E, V, ^h ) (Hamilton-Jacobi principle): in a way 
analogous to the Lagrangian formalism, now using sections of Tq'-V — > B, and the form @h - So 

we look for sections ipo £ T C (B,V) which are critical for the functional H°(^) := / ip*®h > 

Jb 

for every ip £ T C {B,V). Therefore, it can be proved that the critical sections of the Hamilton- 
Jacobi variational principle are the integral curves of a vector field Xh £ 3£(V) (which is called a 
Hamiltonian vector field for the system ( J U E, V, ^h )) satisfying that: 

i(X ho )n ho =0 , i(X ho )n°^0. 

As above, the second equation means that the vector field is TQ-transverse, and we will take 
this condition in the form 

z(X ho )77° = 1 . 

Now, as {J l E, is almost regular and (il^,!] ) is no longer cosymplectic, the above equa- 
tions have no solution in general. In the most favourable cases, Hamiltonian vector fields for 
{J l *E, V , fiho) can exist only in a subset of points of V, and the most interesting situation is when 
this subset is a submanifold of V . 



Remark 5 The construction of a Hamiltonian system associated with a hyper-regular Lagrangian 
system {J l E,Q,£) follows the same pattern, but now V = J U E, V = T*E, and J-C is a diffeomor- 
phism. Then, the Hamiltonian section is now 

h :=F£o TCT 1 

which is a diffeomorphism connecting J U E and J-C{J l E) (observe that it is just the inverse of fi 
restricted to ^FC{J l E)). (See |*§ for details). 



4.3 Hamiltonian constraint algorithm and equivalence with the Lagrangian for- 
malism 

The problem we want to solve is the following: to look for a submanifold P V, transverse to the 
projection Tq-.V — > B, and a vector field X^ £ X(V,P) such that 

1. [i(X ho )n ho ]y = 0, and [i{X ho )rj } y = 1, for every y £ P. 

2. Xh is tangent to P. 

Item 1 is the compatibility condition, and item 2 is the stability or consistency condition for the 
Hamiltonian equations. 



Once again we apply the algorithmic procedure developed in Section 2.2, with the hypothesis 



there assumed. Thus, following this process, we obtain a sequence of constrained submanifolds 

... J b Pl € ...A Pl ^P ^V^J u E. (29) 

For every i, Pi is called the (i + l)th Hamiltonian constraint submanifold, and this procedure is 
called the Hamiltonian constraint algorithm. 
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The only case interesting for us is when the algorithm ends giving a submanifold Pf (with 
dimPy > 0), which is called the final constraint submanifold for the Hamiltonian problem. Then, 
there exists a vector field X^ G X(P, P), tangent to Pf, such that 

[i(x ho )n ho ]\ Pf = o , [iCXboto ]^ = 1 • (3°) 

Now, consider the sequences of dynamical Lagrangian constraint submanifolds (|2~l|) and of 



Hamiltonian constraint submanifolds (^j). Then, we have (see [32]): 



Theorem 5 At every level i = 1, . . . , / of the Lagrangian dynamical constraint algorithm and the 
Hamiltonian constraint algorithm we have: 

1. FCj\{Pi) = Mi (and hence JX;_i(M;) = Pi). 

2. The induced mapping by TC-i-i from Mi to Pi, denoted by TLi'.Mi — » Pi, is a submersion. 

3. For every y~i £ M { , T C^ 1 if Li(yi)) = P£ 1 (J r £o(y i )). Hence 

ker (.FA.)* = ker (JFP*)* = (ker fi £ n X V{lTl) (J 1 E)) Vi . 

4- There is a basis {Q } C C°°(J 1 E) for the set of ith- generation dynamical Lagrangian con- 
straints which is J- 'Co-relat 
tonian constraints; that is, 



straints which is T 'Co-related with a basis } C C°°(P) for the set of ith- generation Hamil- 



Thus, the codimension of Mi (as a submanifold of J^E) coincides with the codimension of Pi 
(as a submanifold ofV). 



( Proof ) We make the proof for the first level (i = 1). For the other cases it follows the same 
pattern. 



1. First we prove that FCo(Mi) C Pi. In fact, for every y G Ml, we have that there exists 
Xy G Ty^E such that 

i(Xy)(n C )y=0 , i(Xy)rjy = l. 

Thus, if TCoiy) = y and Xy = (FCo^y^y), then using FCqQ.^ = and TC^rf = n, we 
obtain that 

i(Xy)(Q ho )y=0 , i(Xy) V l = l 

which implies that y G Pi. 

Conversely, for every y G Pi, the fiber J r £g 1 (y) is entirely contained in Mi. In fact, if 
y G FCq (y) and Xy G T^P is such that i%~ ^h (y) = an d ix- V^iv) = 1) then since TL§ is 
a submersion, there exists Xy G TyJ^E satisfying (,F>Co*)y(^) = Xy and, reasoning as above, 
we conclude that ix s ^c{v) = 0> ^x 9 Vy = 1- Thus, y G Mi and therefore TCq 1 (Pi) Q Mi. 

As a consequence of both results we obtain that TLq 1 (Pi) = Mi, and consequently J 7 £q(Mi) = 
Pi- 
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2. For every y G Pi, as J~Cq is a submersion, then there exists an open neighborhood Uq C V 
with y G Uq, and a differentiable mapping so : C^o ~ * ^-E, such that TCq osq = ldjj . We take 
U\ = Uo (~l Pi, and = so \ui : U\ — > J X E is a differentiable mapping. 

If y' G C/i, then there exists y' G Mi such that J-"£o(y') = so T Co(s\{y')) = .F)Co[ s o(y')] = 
jf = F£ (y')- Therefore, s x (y') G F C^ 1 {F C G {y')) C Mi. That is, si: t/i -» J 1 ^ takes values 
in Mi, and because Mi is supposed to be an embedded submanifold of J^E, then S2 - U\ —* M\ 
is also a differentiable mapping. Moreover, TL\ o s\ — Id^. Then, s\\U\ — » Mi is a local 
section of .T-Xi: Mi — > Pi in y G Pi, and therefore J-C\ is a submersion. 

3. Using the first item we deduce that J r C^ 1 (J T Ci(yi)) = FCq 1 (J r £o(y~i)) and therefore 

ker(jF£j*)y. = ker(jF£o*)j/ 4 - 

, it was proved that 

ker(^£ *) = kerfir n X^V 1 ^)- 

4. We have the following diagram 





V 


T*i 


TJl 


Mi ^ 


Pi 



Consider the basis {Q } of lst-generation Hamiltonian constraints, we must prove that the 
set {FCq^P = cj 1 ^} is a basis of 1 si-generation dynamical Lagrangian constraints. First 
observe that, as J*Q = 0, we obtain that 

= TCUif = itFC^f = iic® 

and thus the functions are dynamical Lagrangian constraints. Furthermore, the number 
of independent lsi-generation Hamiltonian constraints is dim V — dim Pi = r. If m is the 
dimension of the fibers of J- Co, then dim J^E = dim J l *E = dim P+m, and as a consequence 
of the items 1 and 3, dim Mi = dim Pi + m. Therefore, the number of independent lst- 
generation dynamical Lagrangian constraints is 

dim J l E — dim Mi = dim P + m — dim P\ — m = r . 

Notice that the constraints are independent because TC is a submersion and thus J~Cq is 
injective. ■ 



Furthermore, in [32 



In this way we can draw the diagram 



M = J 1 E V = P Q ^ J U E 

In T 31 

Mi ^4 Pi 

Ul T A 



t /-2 f f-2 

I V-1 I Jf-1 

M f-1 Pf-1 



(31) 
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Finally, for the vector fields solutions of the Lagrangian and Hamiltonian dynamical problems 
we have: 



Theorem 6 Let Mf be the final dynamical constraint submanifold for the Lagrangian problem, and 
Pf is the final constraint submanifold for the Hamiltonian counterpart. 

1. Let X c G X^EjMf) such that: 
(a) Xc is tangent to Mf. 



(b) Xc is a solution of the dynamical Lagrangian equations (21). 

(c) Xc is T C-projectable. 

Then TLf*Xc = X^ Q G 3£(J E, Pf) is tangent to Pf and it is a solution of the Hamiltonian 
equations $3(\) . 



2. Conversely, if X^ is a vector field tangent to Pf and it is a solution of the Hamiltonian 
equations ftsdj), then there exists a vector field Xc on Mf such that TCf*Xc = X^ , and it 



is a solution of the dynamical Lagrangian equations (2&). 



( Proof) The results follow using p^), that J-Cf : Mf — > Pf is a surjective submersion and the 
fact that TC* rf = n (see ||). ■ 



In this case, the integral curves of Xc are not, in general, the solutions of the Euler-Lagrange 
equations, since they are not holonomic (Xc is not a SODE, necessarily). The problem of finding 
a solution satisfying the SODE condition, that is, an Euler-Lagrange vector field, is studied in 



Section 4.5 



4.4 Intrinsic characterization of Hamiltonian constraints 

We want to give an intrinsic characterization of the Hamiltonian constraints which define the 



constraint submanifolds Pi, and hence we will apply again the results given in Sections 2.1 and 
^L% . First we take a connection V in Tq : V — > B or, what is equivalent, a fj-transverse vector field 
y G 3i(V), which can be selected such that i(y)rf = 1 holds. 



The connections V in Tq : V — > B, and V in -zr 1 : J l E — > B (given in Section 3.2) can be chosen 



.F£o-related. This means that the vector field 3^ G X(J 1 £') associated with V can be selected in 
such a way that T L^y = y. This is always possible, since T 'Co is a surjective submersion. From 
now on this fact will be assumed. 

In this case, V can be constructed starting from a connection V in tt: E — > B, as in the 
Lagrangian case. 

Observe also that since FC^rf = rj, if i(y)rj = 1, then i(y)r] = 1 too. 
Now, let 7h G Q l (J u E) and u^o G Q 2 (J U E) be the forms defined as 

7h := i(y)tth , Wh : = ^h - r?° A 7 ho . 
Using (^) and after the comment in the first item above, it is obvious that 

FCoUha = uj c , ^oTho = 7£ • (32) 



M. de Leon et al, Singular Lagrangian systems on jet bundles. 



26 



Suppose that y» G Mj (respectively, yi + \ G Mj + i) and that yj = TC^yi) (respectively, yi+i = 
.T-XoG/i+i))- Using (|32]) and the fact that FC^rp) = 77 = (7t l )*(w), we deduce that 

(^£o.) fii (TiM) = T^P< 

(^W^T^M.nT^M.+O = T^ +i pnT 5i+1 p +1 . 

Moreover, from Lemmas ^ and [l?] (see Appendices [A| and |C|) and Theorem ||, we obtain that 

ker(.F£o*)w C (ker ft £ n ker 77)5. C T^M C T^M;, 
ker(^£ *) g! C T^- i+1 M n Ty i+1 M i+1 C T^ +1 Af< D T Si+1 M i+1 . 

Thus, if m is the dimension of the fibers of J-jCq, it follows that 

dim T± Mi = dimT^P + m 
dim(T^ +1 M,nT5 i;+1 M m ) = dim(T^ +i P n T 5i+1 P i+ i) + m. 

Therefore, we conclude that the rank of the distributions T^Mj and T^. i M{ n TMj + i is constant 
if and only if the rank of the distributions T p and Tp. +i Pj n TPj + i is constant. We will assume 
that the distributions T Pj and Tp . x P n TPj + i have constant rank. 

Denote by X" L (Pj) the set of sections of the vector bundle T^Pj — > Pi- Then, as a direct 
application of Proposition [l] and Theorem [3| to the present case, we have: 

Theorem 7 Let (J E, Qh ) be the Hamiltonian system associated with an almost regular La- 
grangian system, and consider the sequence {Pi} of submanifolds ( 2§ ). For every i > 1, let ru_i\, r< 



be the ranks of the distributions T P{-i and Tp.Pj_i n Tp respectively, and consider a local basis 
for X" L (Pj_i) (in a neighbourhood of a point x G Pi) 



1^-1 >--->^'r (i _ij)' a 'r w _ 1) +l»'"'' 2, r j J 



where {Z^ 1 ^}j =1) _ )T . _ x are tangent to Pi. Finally, let X\ M) = + ^ be the general solu- 

tion on the submanifold Pi, where 
parts of the solutions, respectively. 



(i—l) (i— 1) 

tion on the submanifold Pi, where X^ Q and Y h() denote the determined and the undetermined 



1. The submanifold P\, where the Hamiltonian equations are compatible can be defined as the 
zero set of the so-called lsi-generation Hamiltonian constraints, C C°°(V), which are 
characterized as 

£(1) = i(Z^)(rj° - 7h ) , for every £(°> G X ± (V) . 

2. Every submanifold Pj + i in this sequence can be defined (in Pi) as the zero set of the so-called 
(i + l)th-generation Hamiltonian constraints, C C°°(P) ; which are characterized in 
the following equivalent ways: 

(a) For every Z® G X ± (P), 

e (m) =^W-7 h0 ). 

(b) For every j G {1, . . . ,r (i _ 1} }, 

i [l+1) = Uxt l) Mzf~ l) ){^- lh0 )]. 
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(c) For every j E {1, . . . , r(j_ij}, 

=^([4r 1) » i r i) ])^o +L(3>)[z(^- i) ) 7 h ] . 

5. For every k S {r7j_i) + 1, . . . , rj} i/ie stability condition 

o = + ^Mif 1 ^ - 7h )]|p 

determines (partially or totally) (on Pi). 



4.5 The second order differential equation problem 

Assume that ( J X E, tic) is an almost regular Lagrangian system and that Mf is the final constraint 
submanifold (in the Lagrangian setting). Then there exists a vector field xjr on Mf such that 
i(X^)Qc\M f = 0, and i(x£)rj\M, = 1- But will not in general be a SODE on Mf, and thus the 
integral curves of X^ will not satisfy, in general, the Euler-Lagrange equations. In order to solve 
this problem, we will construct a submanifold S of Mf where there exists a unique vector field X^ 
such that X^ is a SODE on S and 

i(X[)n c \ s = , i(Xi)n\ s = l. 

From Theorem ^j, we know that we can choose the vector field X^ on Mf such that it projects via 
TLf onto a vector field on Pf. Then we consider the subset S of Mf defined by 

S = {xeM f /(J(X f c ))(x) = 0}. (33) 

In H^], it was proved that 

J(X f c )(x) E ker^vC/)^ = ker(^£ )*x, for every x E Mf (34) 

and that for every x E Mf, S Pi FCj 1 (FCf^)) = S n -P£q 1 (J r £o(x)) is a point s x f(x). This point 
is characterized by the condition 

i( Sxfc (x)) = (^U(X f c (x)), 
where i : J^E — > TP is the canonical embedding (see Section |3.1| ). 

The above result allows us to introduce a well-defined map s„/ : Py — ► My such that S = 
s x f(Pf) and .PC/ o s x f = Id. In fact, sw : Py — > My is a global section of the submersion 
TCf : Mf — > Py and therefore S is an embedded submanifold of My and the map s x f : P/ — > S 1 is 
a diffeomorphism (for more details, see |32fl). 

Moreover, we have 
Theorem 8 J||/ 

1. There exists a unique vector field X^ = (s x f)*(X^ o ) tangent to S which satisfies the following 
conditions 

i(Xi)n c \ s = , i(Xl)ri\s = 1 , J{X S C )\ S = 0. 
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2. The integral sections of satisfy the Euler- Lagrange equations. 

Next, we will give a local description of the submanifold S and of the vector field X^. 

If m is the dimension of the fibers of the submersion J-Cq : ^E — > V then it is clear that the 

rank of the partial Hessian matrix ( q v q p J is n — m (see (p6|)). We can suppose without loss of 

the generality, that the first n — m rows of this matrix are independent. With this hypothesis and 
using (p6|), we deduce that for every j G {1, . . . , m} there exist local real functions {Wj}i<i;<n-m 
such that {Wj}i<j< m is a local basis of ker(.F>Co)*! where 

W > = e^^ + W iW> f-j €{!,..., m}. (35) 



Now, if 

dt dq % dv l 



then, since X^ is .7-Xj-projectable, it follows that the functions A % are constant on the fibers of 



FCf : M f -» P f . But, as FCl^tfCfix)) = T Cq 1 {T C {x)) for every x G Mf (see Theorem |), we 



obtain that 

W j \ Mf (A i )=0, for j G {1,... ,m} and % G {1, . . . ,n}. (36) 

Furthermore, 
Thus, from ( p4| ) and (|3^), we have that 



J{X f c ) = (A n ~ m+ i - v n - m +i)Wj\ Mr (37) 



Note that the functions 
are independent on Mf. In fact (see (33) and fl36D), 



£S = A n-m+j _ v n-m+j^ I < j < m 



dtf(Wi\ Pf ) = -8, 



I.J- 



Moreover, using ( |33[ ) and (|37|), we conclude that {Cj}i<j<m is a set of local independent con- 
straint functions defining S as a submanifold of Mf, that is, 

S = {i£ Mf/(A n ~ m+j - v n ~ m+ i)(x) = 0, Vj G {1, . . . ,m}}. 

Finally, a direct computation proves that the vector field on S is given by 

Xi = X f c \ s - X^si&WAs = X f c \ s - (X f c \ s (A n - m+ i) - A n - m+ i)W 3 \ s . 

5 Finding Euler-Lagrange vector fields. Dynamical and SODE 
Lagrangian constraints 



In Section 4.5 it was shown that there exists an embedded submanifold S of Mf and a unique 
vector field tangent to S which is a SODE and verifies the Euler-Lagrange equations. However, 
there is another different approach, consisting in imposing the SODE condition on each step of the 
constraint algorithm. Next we will develop this procedure, which is a generalization of the method 
given in |58| for autonomous Lagrangian systems (see also ||). 
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5.1 SODE condition: SODE first order generation constraints 

In Section [3| we have applied the results of Section ^ to obtain a constraint submanifold Mi such that 
there exist vector fields {xfi + Y^ , G ker 0^ n ker n} satisfying the dynamical Lagrangian 
equations (|l]) or, equivalently, 

i{xf ] + yP)u c \ Mi = -1c\m, , i(xf + YP) V \ Ml = 1 . (38) 

However, in general, these solutions on the constrained submanifold M\ do not satisfy the SODE 
condition on M\. Hence, we will look for the points in M\ where such a solution exists. 

Thus, we consider the submanifold 

Si = {x G Mi | 3?P G ker n c n kerr/ such that (J{xf> + Y^ 0) ))(x) = 0} . 

That is, Si is the maximal set of points of Mi where there exist a vector field D = X^ + Y^ G 
3L{J l E; Si) such that 

(a) i(D)uc\ Sl = -Jclsj, i(D)rj\ Sl = 1, and 

(b) D satisfies the SODE condition on Si (J(D)\ Sl = and i(D)f]\ Sl = 1). 

Notice that, if D G X(J X E) is a vector field satisfying (a) and (b), and W G keru;£ n ker^, 
then D + W satisfies (a) and (b) too. In particular, since kerf2£ D ker 77 C keru;£ n ker 77, then 
kerQc n ker7rj C kerw£ n ker7r*. Hence, if W G kerQc H ker-zr^, D + W also satisfies (a) and 
(b). Conversely, if Di and D2 satisfy (a) and (b), then Di — D2 is an element of ker Qc fl ker 7r* C 
ker uj£ fl ker tt\ . 

Now, we are going to describe the submanifold Si inside Mi as the vanishing of a family of 
functions. For this purpose, we will use a particular kind of connections: 

Definition 3 A connection V = rj <g> y in 7? 1 : J 1 E — > B is said to be a second-order connection if 
y is a SODE. 

For second-order connections, we will prove the following results: 

Lemma 2 Let V = n®y be a second order connection in the fibered manifold tt : J l E — > B. Then 
the 2-form ui£ = Qc — V A 1c satisfies 

uJc(J(X),Y) + uJc(X,J(Y)) = 0, 

for every X,Y G X^E). 

( Proof) First, notice that for every X G X^E), i(J(X))n = (as J(X) G ker 71^). Moreover, 
we know that for every 1,7 6 3L{J l E) we have ^(Jpf ), Y) + n c (X, J(Y)) = (see §§), hence 

;(JP0)7£ = iCJWW)"/: = -»(X)t(.7(Wlc = o 

since J{Y) = is a SODE). Then, (r? A j c )(J(X),Y) = 0. Therefore, 

w £ (J(X), Y) = n c (J(X), Y) + (77 A tcXJW, *0 = JPO = JW) • 
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Lemma 3 Let M be the subset ofX(J l E) defined by 

M = {Z G X{J l E) | J(Z) G ker^£,} 

and assume that the connection V = rj®y is of second order (J{y) = 0). Then, given X G X(J 1 E) 
verifying i(X)i(Z)uj£ = 0, for every Z G M., there exists a ^-vertical vector field V such that 
X — V G ker uj£ . Such a vector field V is unique up to an element of ker T . 

( Proof) Recall that ker = ker f^nker tt\. First, notice that if V is a second order connection, 
then ker T = ker Q,£ n ker 7r* = ker uj£ n ker i[\ . We know that ker Q,£ n ker irl C ker uj£ n ker irl . 
But if V G kerw£ n ker 7r*, then there exists a vector field [/ G X^E) such that = V (V is 

7r 1 -vertical). Then 

i{v)n c = i(v)u c + i(v)(7 1 /\ 1 c) = -(i(v)i(y)n c )v 

= -d(J(u)) i(y)n c )v = d(u) i(j(y))n c ) v = o 

(since y is of second order). 

Now, given an arbitrary vector field X, the necessary and sufficient condition for the equation 
i(V)u>c = i(X)u>c to be solved for a ^-vertical vector field V G X(J l E) is that for every U G 
X(J l E) verifying i(U) i(V)uj£ = (with V a ^-vertical vector field), then i(U) i(X)uj£ = 0. Now, 
notice that, if i(U) i{V)uj c = 0, for every V, then i{U) i{J(Y))u c = - i{J{U)) i(Y)u c = 0, for 
every Y G X^E); that is, U G M. Therefore, if X verifies i{X)i{Z)uj c = 0, for Z G A4, then 
i{U)i(X)iOc = and equation z(F)u;£ = i(X)uc has solution. Finally, for every couple V\,V2 of 
such ^-vertical fields, their difference is ^-vertical and belongs to kerw£. ■ 

Lemma 4 Let M' be the subset ofX{J x E) defined by 

M' = {Z G 3L{J X E) \ Z G kervf 1 and J{Z) G kerjFX*} 

and Zei V = n^D^ fre a second order connection. Then for every X G ker Tt\ verifying i(X) i(Z')u>c = 
0, for every Z' G M., there exists a ^-vertical vector field V (unique up to an element o/ker.F£*j 
such that X — V G ker ujc ■ 

( Proof ) First, notice that if V is a 7r 1 -vertical vector field, then there exists a ^-vertical vector 
field Y' such that J{Y') = V: if Y G £( J 1 ^) is such that J(Y) = V and D is an arbitrary SODE, 
then the vector field Y' = Y — (i(Y)n)D is ^-vertical and J{Y') = V. Now, the proof is similar to 
that of the above lemma, simply observing that, if i(U) i(V)uc = 0, for every V ^-vertical, then 
i(U) i(J(Y'))oj c = - i{J{U)) i{Y')uj c = 0, for every ^-vertical vector field Y', that is, U G M' . ■ 

Now, we have the following. 

Proposition 6 Let Y G X(J 1 E) be such that X^ + Y satisfies the SODE condition; that is, 
J(X^ + Y) = and i(X^ + Y)rj = 1, and let V be a second order connection. Then 

Si = {x G Mi | (i(Z) i(Y)LO C )(x) = , \/Z G Ai} , 



( Proof) Consider the set C = {x G Mi | (i(Z) i{Y)uj c ){x) = , VZ G X}. Then the proof can 
be done in three steps: 
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1. C is independent of the chosen vector field Y. 

Take T = Y 1 -Y 2 , where Y\ and Y 2 verify that X ( p + Y { are SODE for every i = 1, 2. Then T 
is a 7Ti-vertical vector field. So there exists a vector field J7 G X(J 1 £') such that T = J{U). 
Now, using Lemma |2] we have 

u c (T, Z) = ujc(J{U), Z) = -uj c (U, J(Z)) = , 

since J(Z) G ker T X* = ker Qc H ker 7T* = ker W£ H ker tt\ . 

2. C C Si. 

Let x € C, then i(Y)u;£)(;r) = for every Z G M. From Lemma ^ we know that at 
each point x G C, there exists a 7r -vertical vector v such that Yx — v G kerw£|aj. Hence 
(xj, -* + y) 5 — v = (X^) s + (Yx — v) is of second order on Mi, and it is also a solution on 
Mi, so x G Si. 

3. Si C C. 

Let x £ Si. Then there exists Y^ G ker$7£ n ker 77 C keriij£ n kern such that + )x 

is a solution of the dynamical equations and satisfies the SODE condition on Si. In general 
J{xf + Y^) x ^ for x £ Si, but a new vector field Y can be constructed such that Y x = 
{?P) X , if x G Si, and [^(X^ + r)] s = 0, for every x G J 1 ^. Then (i(Z)i(Y)u c ){x) = 0, 
for every x G Si (since [i(Y)cj£](:c) = [i(Y^°^)a;/;](x) = 0, for every x G Si), so x G C. ■ 



Corollary 3 With the hypothesis of the above proposition, for every SODE D we have 
Si = {x G Mi I (i(Z) i{xf ] - D)uj c ){x) =0, VZ G M} . 

( Proof ) As y = -X^ — D, the proof is trivial. ■ 

Given Y G ^(J 1 ^), such that xf 1 + Y satisfies the SODE condition, for every Z £ M, the 
function 

</> (1) :=i(Z)i(Y)u c G C 00 ^) 

(or equivalently t^ 1 ) = i(Z)i(X^ — D)ujc, for a SODE L> G %(J l E)), is called a lsi-generation 
SODE Lagrangian constraint. 



5.2 General expressions of first order generation constraints. jF£-projectability 

Using a second order connection V to make the splitting of the form fi£, we can give a unified 
description of both the dynamical and SODE Lagrangian constraints. 

Theorem 9 For any SODE D we have 

Si = {x G J l E I [i(Z){i{D)uj C + lc)]{x) =0, VZ G M} . 

The function i(Z)(i(D)ujc + r yc)) i> s called a lst-generation Euler-Lagrange constraints, and Si is 
called the submanifold of 1 si-generation Euler-Lagrange constraints. 
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( Proof ) Let x G My. If D = X^' + Y for an arbitrary vector field Y, then 

[i(Z)(i(D)u c + y c )](x) = [i{Z)(i{xf)u c + i{Y)u c + lc)](x) = [i(Z) i{Y)uj c ]{x) 
which are the SODE constraints defining S\ as a submanifold of M\. 

Next, we prove that, if Z G X ± (J 1 E), then Z G Ad; that is, J(Z) G ker T = keraj£ Piker 7r^. 
In fact, we know that Z G X" L (J 1 i?) if, and only if, 

i(Z)i(U)u £ + (i(U)r,MZ) v )=0, 

for every U G ?Z(J l E). Then Z G M because 

iCL/HJXZ))^ = -i(J(U))i(Z)u £ = (i(Z) v )(i(J(U)) V ) = , 

Taking this into account, we obtain 

i{Z){i{D)u c + 7£) = + i(^)7£ = -(i(Z) V MD) v ) + i(Z) 7£ = i{Z){ lc - v) , 

which are the dynamical constraints defining Mi as a submanifold of J X E. ■ 

Summarizing, we have arrived at a submanifold S± of J X E and a family of vector fields 

{r c := Xf ] + Yj? ] + V (0) ; V W G kerTC,} (39) 

such that 

[i(r c )n c ]\ Sl = o , [i(r c )rj\\ Sl = i , [J(m\ Sl =o. 

Observe that xf ] and y£ 0) are fixed vector fields, whereas G ker .FiZ* is arbitrary. 
The submanifold Si is characterized as 

S l = {xeJ 1 E\ [i(Z)(i(D)io c + lc)}(x)=0, VZ£M}, 
for any SODE D. The constraints that define Si are of two different kinds: 

• Lagrangian dynamical constraints: C (1) := i(^ (0) )(?? - 7c), £ (0) G X ± (J 1 E). 

• SODE Lagrangian constraints: </>W := i(Z)i(y)w £ , Z G M, Z g" ^(J 1 ^), Y G ^(J 1 ^) 
with x£ 0) + Y a SODE. 

We know that the Lagrangian dynamical constraints can be expressed as .7-X-projectable functions 
(Theorem^). But this is not the situation for the SODE Lagrangian constraints, as will now be 
proved. First, the following results are required: 

Lemma 5 For every V G ker PC*, there exists Z' G M' such that J{Z') = V and Z' is TC- 
projectable. 

( Proof) It is evident that if V G ker ttI is .7-X-projectable, then we can choose an ,F£-projectable 
^-vertical vector field Y' with J(Y') = V (in fact, if Y G 3i(J x E) is such that J(Y) = V and 
D G 3L{J l E) is a SODE. Then Y' = Y - (i(Y)rj)D and the only problem is to take the vertical 
parts of Y and D, which are arbitrary). In particular, for every V G ker we have T£*(V) = 0, 
which implies they are ^"/Z-projectable, hence there exists some Z' G M! which is .F£-projectable 
and J(Z') = V. m 
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Lemma 6 If D G X(J X E) is a SODE, then J(\V,D]) = V , for every ^-vertical vector field V . 



( Proof ) This follows by using the local expressions of the vertical endomorphism J , a SODE 
and a -^-vertical vector field. ■ 

Lemma 7 If D £ X(J l E) is a SODE, then [V,D] G M, for every V G FC*. Moreover, 

[ker FC*, D] + ker vr^ = M' . 

( Proof) The first part is a straighforward consequence of the above lemma. For the second part, 
notice that, for every Z' G M', J(Z') is ^-vertical, hence J(Z') = J([J(Z'),D]) G ker T£*, i.e., 
Z' — [J'(Z'), D] is a vertical vector field. ■ 

Now, we prove the non-projectability of the SODE Lagrangian constraints. 

Proposition 7 If is a 1st -generation SODE Lagrangian constraint, then it cannot be expressed 
as a J- C-projectable function. 

( Proof) If is a 1 si-generation SODE Lagrangian constraint, then according to Theorem |9| it 
can be expressed as 

4® =i{Z){i{D)uj c + lc ) , 
where Z G M, Z G" X ± {J 1 E) and D G X^E) is a SODE. 

First, notice that Z 3L^~{J l E) + ker tt\. If V is a ir l vertical vector field then 3Y G %{J 1 E) 
such that J(Y) = Z and 

i{Z)(i(D)uj c + ic) = i(J{Y)){i(D)uj c + i(y)n c ) = - i(Y) i(J(D))uj c - i(Y)) i(J(y))Q c = 

since D and y are SODEs. Therefore, if Z G X" L (J 1 -E') + ker7r*, then we obtain a Lagrangian 
dynamical constraint. As a corollary, Z G" ker7r*. 

Moreover, Z can be chosen as an element of M' (i.e., Z G M and Z G ker7fj) .7-X-projectable. 
If Z is a SODE, J{Z) = and Z G M But 

i (Z) (i(X>)w £ + 7 £) = i(Z) (i(D)Qc ~ i (D) (t? A 7 £) + 7r) = i (Z) i (D)Q C ~ %{Z)lC + %{P)lc + %{Z)lc) 

= i(z) i(D)n £ + i(D) i(y)n c = i(D) i(y - z)n c = o 

since y — Z is ^-vertical. If Z G" ker ff J and it is not a SODE, then the SODE Lagrangian constraint 
(jP-> can be expressed in an equivalent way as 

=i{Z'){i{D)u c +ic) , 

where Z' = Z — (i(Z)rj)D is ^-vertical; that is, Z' G M' . Finally, from Lemma |B| we know 
that vector fields in M.' can be chosen as .F/I-projectable (notice that .F£-projectable and TL- 
nonprojectable vector fields in M! give rise to the same constraints, since they differ on a vertical 
vector field). 
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Now, suppose that <jft) is .T-X-projectable. Then for every V G we have 

= V(^) = V(i(Z')(i(D)u£ + ic)) = V[i(Z%(D)n c -i(D)( V A^c)+lc)} = V(i(Z') i(D)fl c ) , 
since V[i(Z')(j c - i(D)( V A 7jC ))] = V[{i{Z'){ lc ~ 7/0)] = {Z is vf 1 vertical). 

Therefore, 

o = v(cp^) = v(i(z') i(D)n c ) = -v(i(D) i(z')n £ ) = - i([v, d\) i(z')n c - i(D) uv)(i(z')n c ) 

but since Z' has been chosen to be .F£-projectable and Qc is also .TvX-projectable, then 
L(V)(i{Z')Q c ) = 0. Therefore, 

= V(cf>V) = -i([V,D])i(Z')nc. 

Again using Lemma [7|, we know that for every X' in Ai', X' = [V, D] + W, for some V G ker TL* 
and W G ker nl , so we have 

i(x') i{z')u) C = i(x') i(z')n c - i(x') i(z')(r) a 7£ ) = i(x') i(z')n c = 

= i([V, D}) i(Z')Q c + i{W) i(Z')n c = i([V, D}) i(z')n c = o 

since i(X')n = i(Z')r] = (they are vf ^vertical vector fields) and i(W) i(Z')Q c = i(J(Y)) i{Z')Vt c = 
-i(Y)i(J(Z'))Q c = for some Y G X(J l E) (because J(Z') G kerJFX* C kerfi £ ). Hence we con- 
clude that i(X') i(Z f )ojc = 0, for every X' G M' , and therefore, from Lemma |||, there exists a 
^-vertical vector field V such that Z' — V G kerw/;. That is, Z' G kerw£ + ker7r*. But, for 
^-vertical vector fields, kercj£ = X ± (J 1 E), so Z' G %~ L (J 1 E) +ker7r^, which contradicts the 
hypothesis. ■ 

5.3 Stability condition: new generations of constraints 

In general, none of the fields of the family ( p9[ ) of solutions on S\ is tangent to this submanifold. 
So we must search for the points of S± where a vector field G ker T '£* exists such that 

r c = xf ] + Y^ 0) + V® is tangent to Si. Then we define: 

S 2 = {x G 5i | 3V<® G ker TL* , {xf + Yf + V^) s G T S S X } 
= {x G Si | 3y (0) G ker J= X* such that 

[L(4 0) + ?P + yW)C (1) ](x) = 0, VC (1) , and [L(Af } + Y^ + = 0, V</> (1) } 

where we recall that the lst-generation Euler-Lagrange constraints are given by 

C (i) . = i(z^)( V - lc ) , VZ(°) G X ± (J 1 E) , and 

0« : = i(Z) i{Y)uj c , VZ G M, Z <£ X L { J 1 ^), Y - G £( J 1 ^), x[ 0) + Y SODE . 

But, as the Lagrangian dynamical constraints can be expressed as ,F£-projectable functions (The- 
orem D, and V(°) G ker JF£*, using this in the stability condition for these constraints we have 
that h(V)Q^ = 0, and hence this condition reduces to 

[L(xf + Y^ 0) )C (1) ](x) =0 , for every C (1) • 
Then we have two options: 
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1. [L(YP)C {l) ](x) = 0, but [L(xJ. 0) )C (1) ](x) ^ 0, for every x G Si, and for some C (1) - Then we 
obtain new constraints of the form 

(( 2 ) :=L (xJ. 0) )C (1) G C°°{J l E) . 

2. [L(YP)C^}(x) ^ 0, for every x G Si, and for some C (1) . Then we obtain new constraints of 
the form 

0(2) := L (4 0) +yi 0) )c (1) g c™{j 1 e) . 

All these constraints are called the 2nd-generation Euler-Lagrange constraints, and define the so- 
called submanifold of 2nd- generation Euler-Lagrange constraints S2 Si. Recall that we are 
assuming the hypothesis of Assumption [3|. Then we have: 

Proposition 8 1. The following subsets of Si are the same: 

(a) A:={xeS 1 I [L(x£ 0) )C«](x) = 0, VC« | [L(Y^)(^](x) = 0}. 

(b) B:={xeSi I [ i (ZV)( V - lc )](x) = 0, VZW G X L (Si)}. 

(c) C: the submanifold of Si of zeros of 2nd- generation T C-projectable Lagrangian con- 
straints. 

(This means that the functions {C^} defining this submanifold are the 2nd-generation dy- 
namical Lagrangian constraints^. 

2. The functions (f>^ cannot be expressed as J- C-projectable functions, and they are called 2nd- 
generation SODE Lagrangian constraints. 

( Proof) 
1. B = C: 

In order to prove this we need the following: 

Lemma 8 If jm-M J^E is a submanifold defined by J- C-projectable constraints, and 
js- S M is a submanifold defined in M by non J 7 C-projectable constraints, then X" L (M) = 
X ± (S). 

( Proof ) Recall that 

X ± (M) := {Z G X(J l E) I f M [i(Z)u c ~ (i(Z)n)rj] = 0} , 
X ± (S) := {Z G X(J l E) I f s f M [i(Z)ujc ~ (i(Z)r))r]] = 0} . 

Then it is obvious that X ± (M) C X ± (S). 

Conversely, if Z G X" L (S'), then by definition f M [i(Z)u>£ — (i(Z)rj)rj\ G Q l (M) is a constraint 
1-form for S, which is non .F£M-P r ojectable. Therefore, for every X G X(M, S) being tangent 
to S we have 

3* S i(X)b*M(i(Z)o>C-(i(Z)v)ri)]=0. 

Now take a local basis of X" L (M) consisting of vector fields {X{, Yj}, where {X{\ are tangent 
to S, but {Yj} are not. As M is defined in J X E by .7-X-projectable constraints, then all the 
vector fields tangent to the fibres of TC are tangent to M too. On the other hand, as S is 
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defined in M by non ^^C-projectable constraints, these constraints remove degrees of freedom 
in the fibres of J-C, and hence {Yj} are vector fields tangent to those fibres transverse to S. 
As a consequence 

3* S i(Y J MMZ)uJc - (i{Z)rj)rj)\ = f sJ *M Hjm*^) [i{Z)u c - (i(Z)r,)r,] = 

since i{jM*Yj) i(Z)ojc = 0, because jM*Yj G ker TC* C ker U£, and i(jM*Yj)r] = 0, as rj is 
a 7f 1 -semibasic form and jM*Yj are ^-vertical vector fields. So, we have obtained that, for 
every X G X(M) 

3* S i(X)[jl I (i(Z)uc-(i(Z)r ] )r ] )]=0, 

therefore f M [i(Z)u c - (i(Z)r))r]] = 0, and thus Z G X ± (M). Hence X ± (S) C X ± (M). m 

Taking this into account, we have that B is just the submanifold of 2nd-generation dynamical 
Lagrangian constraints (see item 2. a of Theorem and then the result follows from item 4 
of Theorem |5|. 

CCA: 

As a consequence of Theorem ^, we can take a particular ,F£-projectable solution G 
X{J 1 E). Then, as are lst-generation dynamical Lagrangian constraints, they can also 
be expressed as .F£-projectable functions, and hence the result follows. 

ACC: 

As a consequence of item 2.b of Theorem ||, and bearing in mind that Y^ G ker tic H ker rj, 
we have that every constraint defining B is also a constraint for A, and then A C B = C. 

2. It is immediate from the above item. ■ 

Remark 6 Observe that the expression of the constraints depends only on a particular solution 
Xf ] of the La grangian equation, but the expression of (j)^ involves the vector field YjP\ which 
arises from the SODE condition. This fact justifies the above terminology. 

Furthermore, the stability condition for the lst-generation SODE constraints gives (on S±) 

= [ L (4 0) + ?P + V®)4>M](x) = [L(4 0) + YP)4F> + L(y'°^«](i) (40) 
which is a system of linear equations for TA°) , and we have: 

Lemma 9 The system is compatible at all the points of S\. 

( Proof ) Locally we can take a finite set of independent lst-generation SODE constraints, 
<f>X, ■ ■ ■ , (j>h- As we have said, these constraints remove h degrees of freedom on the fibres of J-C 
Then the matrix of this linear system for V^ ' has maximal rank h, and hence the system is locally 
compatible. However, for every collection of local solutions, a global solution can be constructed 
on S\ using a partition of unity on this manifold. Hence the system is compatible at all the points 
of Si. ■ 

Thus from this system we can determine (total or partially) the vector field V^°\ The stability 
of SODE constraints does not give new constraints but removes degrees of freedom in the vector 
fields solution. The solutions can be written in the form 



r c :=xW+YP + vV + V<» 
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where is a solution of @, and G ker TL* is any solution of the system h(V^)(cp^) = 
(on Si), and this contains all the gauge freedom. The solution i~£ is a SODE and is tangent to S\ 
at the points of S2, but in general, it is not tangent to S^. 

This discussion has been carried out for the submanifold Si, but it can be extended recursively 
until no new constraints appear. In every step, we have a submanifold Sj (i > 1), the so-called 
submanifold of ith-generation Euler-Lagrange constraints, which is defined by constraints of two 
kinds: 

• {C }j which are the ith-generation dynamical Lagrangian constraints, and can be expressed 
as the (only) .F£-projectable constraints. 

• }j which are the ith-generation SODE Lagrangian constraints, and are not .7-X-projectable 
functions. 

Now, we must take the corresponding vector field solutions (on Si) 

r c = xf + + 

where G ker TL* is a fixed vector field, and G ker T X* is undetermined. Therefore, 

the tangency condition for these ith-generation Euler-Lagrange constraints leads to similar results 
as those obtained for Si. 

In this way, we obtain a sequence of constrained submanifolds 

..M s/C.^sAjMi^J^E (41) 

and this procedure will be called the Euler-Lagrange constraint algorithm. As will be seen in the 
next subsection, if the final dynamical constraint submanifold Pt exists, this algorithm ends by 
giving a submanifold Sf Pf which is called the final Euler-Lagrange constraint submanifold. In 
such a case, there exists a vector field Xc G %(J l E, Sf) tangent to Sf such that 

[i(r c )n c ]\ S f = , [i(r c )n]\ Sf = 1 , [J(r c )]\ Sf = o . 

As a summary of all the results given in this section, we have proved the following: 

Theorem 10 Let (^E^c) be an almost regular Lagrangian system. Consider the sequence of 
submanifolds and assume that, for every i > 1, the distributions T^Si and T^. i iS , j fl TSj+i 
have constant rank. 

1. The submanifold Si, where there exist Euler-Lagrange vector fields (solutions of the dynamical 
Lagrangian equations satisfying the SODE condition) can be defined (on J l E) as the zero set 
of the lsf-generation Euler-Lagrange constraints, which are characterized as 

i(Z)(i(D)uJc + lc) , for every Z G M, and D G X^E) a SODE . 

These constraints are of two kinds: 

(a) The lsf-generation dynamical Lagrangian constraints {C^} C C°°(J 1 E), which are 
characterized as 

C {1) = i(Z (0) )(n - lc) , for every G X ± {J 1 E). 
All of them can be expressed as FC-projectable functions. 
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(b) The lst-generation SODE Lagrangian constraints {4>^} C C co {J l E), which are char- 
acterized as 

0« := i(Z)i(Y)uj£ G C°°{J l E) , for every Z £M 

or equivalently as 

<j>Q) = i(Z) i{xf ] - D)lo c , for every Z eM 

where Y G 3i(J l E), such that xf ] + Y is a SODE, and D G X^E) is a SODE. 
None of them can be expressed as a T ' C-projectable function. 

2. For every i > 1, the Euler- Lagrange vector fields solution on the submanifold Si can be written 
as 

r c -.= 4 0) + yP + + 

where is a particular solution of the dynamical Lagrangian equations, Y^ G ker uj£ n 

ker rj is a fixed vector field such that J{X^ + Y^ )\s t = 0, V® G ker T X* is also a fixed 
vector field (with = 0), and V® G ker J-C* denotes the undetermined part of the solution. 

3. Every submanifold Si+\ (i > 1) in this sequence can be defined (in Si) as the zero set of the 
(i + l)th-generation Euler-Lagrange constraints, which are of two kinds: 

(a) The (i + l)th-generation dynamical Lagrangian constraints, {C (i+1) } C C°°(J l E), which 
are obtained by making 

C (m) :=L(4 0) )C W 

for every C {i) such that [L(>^ 0) )C (i) ] k = 0. 

All of them can be expressed as J 7 C-projectable functions. 

(b) The (i + l)th-generation SODE Lagrangian constraints, C C°°{J l E), which are 
obtained by making 

^+i) :=L (x( 0) + ^ 0) )c« 

for every C (i) such that [LC^K^Hs, = °- 

None of them can be expressed as a J 7 C-projectable function. 

4- For every ith- generation SODE Lagrangian constraint (jfo (i > 1), the stability condition 

l(xf + ?P + V^-V + \ Si = o 

determines (partially or totally) the undetermined vector field V^~^ (on Si). 



As an evident consequence of this, we have: 



Corollary 4 For every i (1 < i < / ), Si is a submanifold of Mi, and J-C(Si) = T 'C(Mi) = Pi. 



Remark 7 For autonomous almost regular mechanical systems, there is another relation between 
the Euler-Lagrange constraints and the Hamiltonian constraints, which is established using the so- 
called time-evolution operator Q, ||, [22], [23]. The generalization of this relation to the present 
case is in progress. 
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5.4 Properties of Euler-Lagrange vector fields 



Summarizing, in Sections 4.5 and 5.3 we have given two different procedures for obtaining subman- 
ifolds of Mf <^-> J l E where Euler-Lagrange vector fields exist (satisfying the tangency condition). 
Concerning these procedures, the following must be pointed out: 



• The submanifold S in Section |4.5| is constructed from a previously chosen (.FjC-projectable) 
vector field xjr G X(Mj). Hence, in general, there is not only one, but a family of subman- 
ifolds {S} which are diffeomorphic to Pf J 1 * E (and TC{S) = Pf). Moreover, for every 
submanifold S of the family, there exists a unique Euler-Lagrange vector field on S (tangent 
to S) which, in addition, is .F£/-projectable on the points of S. 

• The SODE constraints defining the submanifold Sf as a submanifold of Mf are not TL- 
projectable, and hence they remove degrees of freedom on the fibers of the submersion TLf. 
As a consequence we have that J 7 C(Sf) = FC(Mf) = Pf. Therefore, the submanifolds of the 
family {S} are embedded submanifolds of Sf. Furthermore, the Euler-Lagrange vector fields 
on Sf are not unique. 

Of course, in the particular case where Sf is diffeomorphic to Pf, then the family {S} is made 
up of a unique submanifold which is just S = Sf, and the corresponding Euler-Lagrange vector 
field is unique and .7-Xj-projectable on the points of Sf. 

• Observe also that if dim Sf > dim S, then the Euler-Lagrange vector field solution on Sf is 
no longer ^"£j-projectable on the points of Sf. 

In fact, let x\ = (ti,qi,v\), and X2 = (^2,92,^2) be two different points in Sf, but in the 
same fibre of Tdf. Then, t\ = £2 , and q\ = q2, but v\ ^ V2- Now if D G X^E) is an 
Euler-Lagrange vector field on Sf we have that 



D~ = — 









xt Jl d v p 





d 


^ p 9 








22 ~ dt 


x 2 +V2 dqP 


+ f — 

x 2 + h dvP 


%2 



and FCi^xi) = TCix^), but J 7 C it (Dx 1 ) 7^ FC*{Dx 2 ), and the result follows. 



6 Dirac brackets and time-dependent constrained Hamiltonian 

systems 

In this Section, we introduce the Dirac bracket as a local Poisson bracket on the restricted momen- 
tum dual bundle J U E associated with an almost regular Lagrangian system. Then we prove that 
if g is an observable (that is, g is a C°°(M, M) -differ entiable real function on Pf), the time evolution 
of g consists essentially of two terms: the restriction to Pf of the Dirac bracket of G and a suitable 
Hamiltonian and the restriction to Pf of the derivative of G with respect to a fixed vector field. 
Here G is an arbitrary extension of g to J l *E. 

In order to introduce the Dirac bracket, we use a suitable cosymplectic structure on J l *E and 
consider a special basis of constraints for the submanifold Pf. In fact, we classify the constraints 
into first and second class and show that the second class constraints are Casimir functions for the 
Dirac bracket. 
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6.1 Cosymplectic structures on the restricted momentum dual bundle associ- 
ated with an almost regular Lagrangian system 

Assume that (J 1 -^, tic) 1S an almost regular Lagrangian system. 



We will use the notation of Sections 4.1 and 4.2. The following commutative diagram illustrates 
the situation 




Now, let V be a connection in ir : E — > B or, equivalently, a vector field y E on E such that 
TT*(m)(y E ) = 1. 

For every point x G E, we have a splitting of the tangent space T X E 

T x E = H x (V)®V x (tt), 

where V(ir) is the 7r-vertical subbundle and H(V) is the horizontal subbundle associated with 
V. Denote by Hor^ : T X E — > H X (S7) and by Ver^ : T X E — > V x (tt) the horizontal and vertical 
projectors. Then, 



Hor^(X) = (ir*(m) x (X))y E (x), Ver^(X) =X- Hor^{X). 



(42) 



Thus, if *(Ver v ) : T*E — > T*E is the adjoint homomorphism of the homomorphism of vector 
bundles Ver v : TE -> V(tt) C TE, then it follows that 

\Ver y ) x {a) = a - {%y E ){a)-K*{w) x , (43) 

for x G E and a G T*E, where i y E '■ T*E — > IR is the function defined by 

(i y E )(P) = 0(y E (x)), for G TIE. (44) 

From (|43j), ([44|) and the definition of the Liouville 1-form of T*E, we deduce that 



t(yer v )*(e) = Q-{iy E )(a l Y{m) 
*(Uer v )*(0) = n + d(iy E ) Aia 1 )*^). 



(45) 



Using ( |43] ) it is also easy to prove that *(Uer v ) induces a smooth map *(Ver v ) : J 1 *^ 
T*E" in such a way that 



T*£ 
A T*E 



*(Fer v ) o ^ = *(y e r v ), /i o *(Uer v ) = Id. 



(46) 
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As a consequence, t (Ver^) is a global section of the submersion \x : TE — > J l *E and, in particular, 
= t (Ver^)(J l * E) is an embedded submanifold of T*E of codimension 1. Moreover, the map 
t iVer^) : J l *E — > S E is a diffeomorphism. Note that 

S% = {a€T*E/(iy E )(a) = 0}. (47) 
Now, we introduce the 1-form V and the 2- form cD v on J l *E given by 

@v = t(^v)*(e) : oj w = *(F^v)*(n) = -de v . (48) 

The 1-form O v (resp. the 2-form u) v ) is called the Liouville 1-form (resp. 2-form ) of J l *E 
associated with the connection V. 



Denote by fj the 1-form on J E given by fj = (t )*(tu) (see Section 4.1). 



Theorem 11 1. The couple (uo^,rj) is a cosymplectic structure on J l *E. 
2. If Oy is the Reeb vector field of(oj^,fj) we have that 

t (Ve~y)^ V )=X t y E \ s v, 

where X^y E is the Hamiltonian vector field on the symplectic manifold (T*E,Q) associated 
with the function iy E ■ T*E — > R. 

( Proof) If i\ : S E — > T*E is the canonical inclusion, then since t (Ver^ 7 ) ((a 1 )*^)) = fj, we 
must prove that the couple ((i^)*(f2), (i^)* ((cf 1 )* (w))) is a cosymplectic structure on S E with 
Reeb vector field X i y E \ s v. 

It is clear that (respectively, (i E )* ((a 1 )* (w))) is a closed 2-form (respectively, 1- 

form) on S E and that the restriction of Xiy E to S E is tangent to S E (see (|4^)). Moreover, 
from ( (44| ) and since TT*(w)(y E ) = 1, it follows that ((a 1 )*(w))(X i y E ) = 1, which implies that 
((i%)*((aiy(m)))(X iyE \ sI ) = l. 

Furthermore, using (^), we obtain that ker((i^)*(0)) = (X i y E \ s v). 

This ends the proof of the result. ■ 



Remark 8 From (|5|), (|6|) and (|48|), we deduce that 

m *(gV) = @ -(iy E )((a'y(w)), 

//*(w v ) = O + d(i y E ) A (cr 1 )*^). ^ 

Next we study the relation between the cosymplectic structure (ui^ ,fj) and the cosymplectic struc- 
ture on J U E defined by another connection V' on it : E —* B. 

Suppose that y' E is the vector field on E associated with V. Then, ir*(w)(y' E ) = 1 and 

v = y' E - y E 
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is a vertical vector field with respect to tt : E —* B. This implies that the function % V : T*E —* R 
given by 

(i V")(a) = a(V(a?)), for a G T*.£ 

induces a smooth function f V : J 1 *^ — > R in such a way that 

Al *(iF)=iF. (50) 

Moreover, we obtain that 

Theorem 12 (i) If (cj v ' ,77) is i/ie cosymplectic structure on J l *E defined by the connection V' 
f/ien 

cD v ' = w v + A fj. 

(ii) If f : J^*E — > M is a reaZ C°° -function and X^ (respectively, X^" ) is the Hamiltonian 
vector field of f with respect to the cosymplectic structure (u^,fj) (respectively, ,fj)) then 

Ay -Ay . 

(m) If { , } v (respectively, { , } v 'j is i/ie Poisson bracket associated with the cosymplectic 
structure (vF ,fj) (respectively, (w v ,77) J i/ien 

{ , r' = { , p- 

( Proof) Using @ and it follows that 

^*(w v ') = //(£ V ) + d{iV) A fj = M *(£ V + d{iV)) A 77. 
Therefore, since u* is injective, we conclude that 

£ v ' = 5 V + A fj. 

This proves (i). 

(ii) and (777) follow from (7) and Proposition 10 (see Appendix^). ■ 



Let us now recall the definition of the Hamiltonian density and of the Hamiltonian function 
associated with the Lagrangian system, the 1-form w and a connection V on tt : E — > B (see 



I 13. E 



Using the connection V, we can introduce the Hamilton-Cartan 1-form on V associated 
with the Lagrangian system, the 1-form w and the connection V. This 1-form is given by 

(elUX) = h (x){Ver^ x ({rlU(X))), (51) 

for x € V and X € TxV. The 1-form @7 allows us to introduce, in a natural way, the Hamilton- 



ho 

Cartan 2-form on V which is defined by fi^ = —dO^ Q . 



If X is vertical with respect to the projection Tq : V — » B, a direct computation, using fl42|), 
(^p and the fact that @h = (jo ho)*(&), shows that 

Therefore, there exists /ig 7 £ C°°(V) such that 

©^o " 0^o = ^((fb 1 )*^)) = W). (52) 

Pip 7 = /iQ 7 ( r /°) an d ^q 7 are cane d the Hamiltonian density and the Hamiltonian function associated 
with the Lagrangian system, the connection V and the 1-form w (see || |~J, |l^|). 
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Remark 9 If y' E is the vector field on E associated with another connection V' on tt : E —* B 
and V = y' E — yE is the corresponding vertical vector field with respect to it : E — > B, then using 
( PH ) and (|52|), and the fact that 0/^ = (jo o /i o )*(0), we deduce that 

/^' = ^-iV| P) (53) 

where iV : J^*E —> M is the real function on J U E induced by the vector field V. This implies that 

< = ©I, - Gv)W> < = nZ + d(iV)\p a r? . 



Next, we will introduce an extension hy : J U E -> T*E of the diffeomorphism h : V -> P. 
For this purpose, we consider an arbitrary extension /i v : J 1 *^ — > R of the Hamiltonian function 
/i^ 7 : 7> -> R. Then, we define the map H v : T*E -» T*£ given by 

tf v ( a ) = a _ {h w (n(a)) + (i ^(a))^*^)^^^). (54) 

It is clear that cr 1 o = cr 1 . Furthermore, we have 



Lemma 10 The following relations hold 

(j,oH v = ll (55) 

(i? v )*(9) = e-^o/z + i^)^ 1 )*^) (56) 

(H v y(n) = n + (p*(dh v ) + d(iy E )) A(a x )*H (57) 

H^oFCq = TCq (58) 



( Proof ) From (]5J) it follows (|55| 



Furthermore, using (^4|), the definition of and the fact that cr 1 o = a 1 , we deduce that 
( |56l) holds. Therefore, 

(# v )*(tt) = -d((i/ v )*(e)) = o + (//(<ta v ) + d(iy E )) a (a 1 )*^). 



Finally, we will prove that (^). Let y be a point of J 1 ^. Assume that 7r 1 (y) = ^{TC^y)) = cc 

*.y( 



and that y E {y) is a tangent vector to J l E at y such that (vr 1 ) > , I ^(3^£;(y)) = 3^£i( x )- 



From (H), ©, (HD and (H), we obtain that 

/ l v (M^c (y))) + ^A)(y)(^(x)) = ^(^£ (y))) + e^(y)(^(y)) 



ti*_{TC Q {y)) + e fto (^£ (y))((^ )^(3> E (y))) 
/ lo (^o(y))(^erJ(^(x))) = 0. 



eZ(FCo(m(F£o)*y(yE(m m] 



As a consequence, using (54) and (|59|), we conclude that H^ 7 (!FCo(y)) = !FCo(y). 



Suppose that /u(a) = Lt(a'), for a, a' € T*i£. Then, i? v (a) = H^ 7 (a') (see (|54])). Thus, there 
exists a mapping A v : J 1 *^ — > T*.E such that 

li v o|] = ff v (60) 

Moreover, from (^), we have that (jiio/i 7 )o^ = ^u, which implies that 

fj,oh v =Id. (61) 
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Since \x is submersion, we deduce that /i v is a smooth mapping, and therefore /i v is a global section 
of n : T*E -> J l *E. 

Furthermore, using (pi), (ED and the fact that fj,(P) = V, it follows that h? ' (V) = V. Further- 
more, using (|8|), ( |60|) and fl61|), we obtain that 

Ho o /i v |-p = Id, /i V |-p o i_i = Id, 

that is, hy\-p = Hq 1 = ^0- We now introduce the 1-form G^v and the 2-form Vt h \j on J l *E given 
by 

e ftV = (/ i v )*(e), n h v = (h^)*(n) = -de h v. 

It is clear that 

f (@ hV ) = @ h0 , f (n h v) = n ho . (62) 

From ( |57| ) and (|60|), we deduce that 

/x*(fi h v) = + (H*(dh v ) + d(j A {o l )*{w) 
and using (B^), it follows that 

H*(tt h v) = //(cD v ) + H*(dh V ) A (a 1 )*^) = ^{u 7 + dh? A 
Since u* is injective, this implies that 

n hV = w v + d/i v A fj. (63) 
Thus (see Proposition 10 in Appendix |C]), we have the following 



Theorem 13 1. The couple (Q h v,fj) is a cosymplectic structure on J^*E with Reeb vector field 
lZ h v given by 

TZ h v = 7£ v + A^v = -E/jV , 

where X^, ( respectively, E^ ) is the Hamiltonian vector field ( respectively, the evolution vector 
field) of h with regard to the cosymplectic structure (uo^,fj). 

2. If f ■ J l *E —* R is a real C°° - differ entiable function and X^ (respectively, X^ ) is the 
Hamiltonian vector field of f with regard to the cosymplectic structure (Q h v,fj) (respectively, 
(w v ,r/) then xf = Xj . 

3. If { , } v (respectively, { , ) is the Poisson bracket associated with the cosymplectic 
structure (w v ,r/) (respectively, (Q h v,fj)), then 

{,P = {,} hV . 



Remark 10 If (J l E,VL c ) 

is a hyperregular Lagrangian system then V = J^E, jq = Id and 
Qh = ^^v. Thus, (Qh ,r/°) is a cosymplectic structure on J U E and TZ h v = E^ is the unique 
solution of the Hamilton equations 

i(K h v)n ho =0, i{K h v)if = l. 
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Next we will write the local expressions of some of the geometric structures introduced above. 
Assume that (t,q 3 ,pj) and (t,q 3 ,p,Pj) are natural coordinates on J^E and T*E, respectively. If 
the local expression of 3^ is 



then 



*(yer v )(i,<z\p, Pi ) = t (VerV)(t,q i , Pi ) = (t,q i ,-yi(t,q i )p j ,p i ) 
G v = pidq* - {y^qi^dt 

( dy i 



Cj w = dq % A — dt A I 3^ (t, g J )dpi + P^~§^ d< l 



dt dq 1 dq 3 dpj 

H^(t,q\p, Pl ) = (^q^-hV^q^p^-y^qi^Pi) 

h V (t,q\Pi) = M,-h*M, Pj ) ( 64 ) 
S hV = pM + ^^q^p^-y^q^dt 

( dy i \ 

n h v = dq' A d Pi - dt A ( y% (t, q 3 )dpi + Pi-^-dq 3 \ + dh v A dt 
dh v d dh v d 



X 



x 



hV dpi dq 1 dq 1 dpi 

d i, dh v \ d /dh v dyi\ d 



j 



6.2 First and second class constraints and the solutions of the Hamiltonian 
dynamics 



Assume that (J 1 E,£lc) is an almost regular Lagrangian system as in Section 6.1 



If V is a connection on it : E — > B, we will denote by (o) v ,fy) the cosymplectic structure on 
J l *E introduced in Section |6.1| (see Theorem |ll]), by A v the Poisson 2- vector associated with such 
a structure, and by { , } v the corresponding Poisson bracket (see Appendix |B|). 

Let Po (resp. Pf) be the first (resp. final) Hamiltonian dynamical constraint submanifold. If 
Po (resp. pf) is the dimension of Po (resp. Pf), then since the pull-back of the 1-form fj to Po 
(resp. Pf) is not zero at every point of Pf, it follows that the distribution # Av ((TPo)°) (resp. 
#^v((TP/) )) has constant rank 2n + 1 — po (resp. 2n + 1 — pf) along Pf. We will suppose that 
the distributions TP n # Av ((TP )°), TP f n # A v ((T P/ P ) ) and TP f n # A v((TP/)°) also have 
constant rank ko,kof and kf, respectively, along Pf. 



Remark 11 If V' is another connection on ir : E — > B, then from Theorem 12 we deduce that 
A v = A v '. Thus, 

TP n# Av ,((TP ) ) = TP n# Av ((TP ) ) 
TP/n# A v'((T P/ Po) ) = TP / n# A v((T P/ Po)°) 

TP/n# A v((TP/) ) = TP / n# Av ((TP / )°). 

With the above hypotheses, we will show that it is possible to choose a suitable basis of con- 
straints for the submanifold Pf. We will proceed in three steps. 

First step: Assume that {Ci°^}i=i,...,2n+i-p is a set of local independent constraint functions 
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defining Pq as a submanifold of J U E. Then the rank of the matrix 

U?i 'V / |P / Jl<i,i'<2n+l-p 

is Zq = 2n + 1 — po — ko- We can suppose, without loss of generality, that the Iq first rows of this 
matrix are independent. With this hypothesis, it is easy to prove that for every j, 1 < j < ko, there 
exist local real functions {/j}i<i<; on J X *E such that the matrix 

(rA°) t(°)-i"V| \ 

USi >?i' I |P/Jl<i,i'<fo 

is regular. In addition, }i<j<k is a local basis of the distribution TPo n #^ v ((TPo)°) 

along Pj, where 

«f=C-/j<f. for !<,'<*„. 

Thus, if ^ = for 1 < i < Iq, it is clear that {Ci°\if^} i<i<i a is a set of local independent 

i<j<k 

constraint functions defining Pq as a submanifold of J l *E. 

Second step: Assume that , ^ , ^ } i<i<i is a set of local independent constraint 

1 < j < k 

1 < r < p a — pf 

functions defining Pf as a submanifold of J l *E. Then the rank of the matrix 

I<r<p -P/ 

is ko — kof. We can suppose, without loss of generality, that the ko — kof first rows of this matrix 
are independent. With this hypothesis, it is easy to prove that for every t, 1 < t < kof, there exist 
local real functions {ff[ }i<j<k -k of , such that 

{ x go)\p f }i<t<ko f 
is a local basis of the distribution TPf n #^ v ((Tp f Po)°), where 

ft - Sko-k of +t - 9tSj > for 1 < t < fc 0/ . 

It is clear that {CfMfM}^,^}, with 1 < i < Z , 1 < j < k - k of , 1 < * < k of and 1 < r < 
Po — Pf, is a set of local independent constraint functions defining Pf as a submanifold of J 1 *P. 

Third step: We consider the (Iq + (ko — kof) + (po ~Pf)) x Co + (ko — &o/) + (Po ~Pf)) matrix 

A direct computation shows that the Iq + (ko — kof) first rows are independent and that the rank 
of this matrix is Iq + (ko — kof) + Sf, with Sf = (po — pf) — (kf — kof)- We can suppose, without 
loss of generality, that the Iq + (ko — kof) + Sf first rows are independent. Then, it is easy to prove 
that for every u, 1 < u < kf — kof, there exist local real functions {h r u }i< r < Sf such that 

\p f , x~y f) \p f } i < t < k 0J 

St Su 1 < u < k f - k f 
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is a local basis of the distribution TP/ n # Av ((TP/)°), where 

Kf) = _ h rdf) 

In conclusion, if & f) = , for 1 < r < s f , we have that {g 0) , if , 4 (0) > Cr /} , ft/* } with 1 < i < Z , 
1 < i < ko — k f,l < t < k j, 1 < r < Sf and 1 < u < kf — koj, is a set of local independent 
constraint functions defining Pf as a sub manifold of J l *E. Moreover, along Pf 

T-Po n #Jvv((TPo)°) = ( x Jo)\p f , x J )\Pf) i<j<k -k of 

?j St 1 < t < k f 

TP f n#~ AV ((T Pf P )°) = (X^Ip^^ (65) 
TP / n# Av ((TP / )°) = (XJ^Ip^X^Ip,) 

St Si" 1 < u < kf — k f 

In addition, if we denote by 

(66) 

1 < r < s j 

then we deduce that the matrix 

(C a p = {X a ,X/3} V ' \p f )l<a,p<l +(ko-k of )+s f 

is regular. 

Following Dirac's terminology, {X a }i< a <i +(k -k of )+s f i < t < k of are said to 

1 < u < k j — kQj 

be second-class and first-class constraints for the submanifold Pf, respectively. 

We now consider the completely integrable distribution D whose annihilator D° is generated 

by 

{dX a }l<a<l +(k -k Qf )+Sf 

Since the matrix {C a p) is regular, it is follows that there exists an open neighbourhood U of Pf in 
J U E such that TV J U E = D® # Av (D°). 

Thus, we have the corresponding projectors 

P : T v J U E -► D, <Q> : T V J U E -» # Av (P>°). 

A direct computation shows that 

Q = C a/3 Xj gxity, P = Id - C a/3 Xj a <g> dAjj, (67) 

where (C a/3 ) is the inverse of matrix (C Q( g). 

Remark 12 If X € X(U), we have that 

cD v (P(^ v ),P(X)) = -w v (Q(ft v ),P(X)) = -C a ^1l v {Xp){i{Xj)0j v ){¥{X)) 

= (c^vy{x p )vy{x a mnx)). 

Thus, using that C af3 = -& a , it follows that 

cD V (P(^ V ), P(X)) = 0, r?(P(^ v )) = 1. 
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Therefore, if TPf n ((TP/) ) = {0} and jf : Pf — > J 1 *^ is the canonical inclusion then 
P(7^ V )| P/ is tangent to P f and 

*(P(£ V )| P/ )#(£ V ) = 0, i(F(>R?)\ Pf )f f (r}) = I- 

As a consequence, from Lemmas |TJ] and [16] (see Appendices [5| and |B|), we deduce that the couple 
(yj-(d) v ), J*f(fj)) is a cosymplectic structure on Py, and that P(7£ V )| P/ is the Reeb vector field of 
this structure. 



Remark 13 If V' is another connection on tt : E — > B, and we consider the same initial set of 
constraints {Q }i=i,...,2n+i-p defining Po (as a submanifold of J U E), then following the above 
process we obtain the same final set of constraints {Q ' , £j , ^ , £^ , £u } defining the submanifold 
Pf (see Theorem |l~2| ). Therefore, we also obtain the same splitting of TjjJ^E and the same 
projectors 

P : T v J U E - £>, <Q> : T V J U E - # Av , (P>°) = # Av (£>°). 

Next, suppose that /i v : J l *E — ► K is an arbitrary extension of the Hamiltonian function 
/ig 7 : Po ^ M and that = 7£ v + is the evolution vector field of h7 with respect to the 

cosymplectic structure (a) v ,f/) (see Section Bj]). 



We will prove that the restriction to Pf of F(E^ 7 ) is tangent to Pf, and that such a restriction 
is a solution of the constrained Hamilton equations. For this purpose, we will use the following 
result. 

Lemma 11 If xf is a point of Pf, then 

{X G T Xf P /i(X)(n ho (x f ))\ TxfPo = 0, i(X)(r,\x f )) = 0} = T Xf P Q n # Kv (x f )((T Xf P )°). 

( Proof) Suppose that X^ o is a solution of the constrained Hamilton equations on the submanifold 
Pf, that is, X[ Q G X(P/) and 

(^4)^)^=0, (i(X f Jr,°)\ Pf = 1. 

We consider the cosymplectic structure (iD v (x/), rp(xf)) on the vector space T Xf J l *E, where 

cD v (x / ) = n h0 (x f ) - (i(Yl)n h0 )(x f ) Ary ^). 

The Reeb vector of this structure is Y^(xf). Since Y^(:ej) G T Z/ .Po and 

a) v (x / )(F,Z) = ^(x / )(y,Z), 

for T Xf Po, it follows that 

{X G T^Po/WXXO^x/)))^ = 0, i(^)(r7 ^/)) = 0} 
= {X G T^PoMAO^O^))!^ = 0, ipO^s,)) = 0} 
= T Z/ P fl #av( X/ )((T I/ Po)°), 

where A^(xf) is the 2-vector on T Xf J u E associated with the cosymplectic structure (ui)^(xf), 
fj{xf)) (see Appendix H). Therefore, using Proposition || (see Appendix || ) and the fact that 

Cj V {x f ) = ^(xf) + (dh V (x f ) - (i(Y/jn h0 )(x f )) A fj(x f ), 
we deduce the result. ■ 
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Theorem 14 The restriction to the submanifold Pj of the vector field F(E^) is tangent to Pf 
and, in addition, this restriction is a solution of the constrained Hamilton equations; that is, 



[i(¥(E^))Q ho ]\ Pf = 0, [i(F(E^)m Pf = l. (68) 

( Proof ) The definition of the distribution D implies that 

F(E^)(X a ) = 0, for every a. (69) 

Furthermore, we have that 

F(E% ) = E$- A'XTo) - ^Xj 0) - SXln . (70) 

?i ?j Sr 



(71) 



As a consequence, from ( |65| ) and (frCp, we obtain 

P(£$)|p,(g 0) ) = ^| P/ (|f)), 1 < t < k f 
F(E~£)\p f (& f) ) = e J\pA8P)> !<«<*/- ko f . 

Now, let be a solution of the constrained Hamilton equations on the submanifold Pf . 

Using dH), (§3^) and the fact that 

n h0 (x f h0 ,xp\p f = n h0 (x[ () ,xp\ Pf = n hQ (x[ o ,xp\ Pf = o 

we deduce that 

^vMff ) = E^\ Pf (^) = E,*\ Pf m = (72) 
for 1 < j < ko — kof, 1 < £ < ^o/ an d l<w</c/ — /cq/- 



Therefore, from fl69|), (|7l|) and (|72|) , we conclude that the restriction of F(E^y) to is tangent 
to Pf. 



Next, we will prove that ( pq) holds. Using (p2j), (p3|) and (|T2|), we have that 

?(0) t(0) lf V 



W^m)^)!^ = °j for ever y i> i < i < - fco/- (73) 



Furthermore, since the matrix ({^^HIp,) 

l<i,i'</ is regular, one can find local real functions 
{^*}i<i<Jo on J^E such that the restriction to Pf of the vector field 



y0 ZP~V \i V~V 

is tangent to Pq- Moreover, it follows that (see (62) and Theorem 

(i(x° ho )n ho )\ Pf = (#£$($V)Ip,- 

But as ^ |p is tangent to -Fbi we have 

= (VAt4 0) ))| P/ = (X j E^0 ] ))\p f ~ (A J P'{ef ,Cf 0) }" v )|p / = (V^(^ 0) ))|p r 
Thus, we deduce that 

(i(X° hQ )n h0 )\ Pf = 0, (»«)>) )|p y = 1. (74) 
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This implies that 

(»(*4 - x° ho )n ho )\ Pf = o, (i(x[ o - X° h0 )n )\ Pf = 

and consequently (see ( |65| ) and Lemma [ll]) 

X[ Q = (E h l - A*Xj } - tfX^ 0) - R k Xj 0) )\ Pf , (75) 

where /2 J and K fc are local real functions on J^*E. 
that the restriction to P/ of the vector field 



Now, from (|70|), (|75|), and since (X^)^ and (F(E^ 7 ))\p f — x[ are tangent to Pf, we obtain 



z = (\ l - V)x7 0) + - - ^^J/) 



is also tangent to Pf. Hence, it follows that Z(xf) € D(xj) n #^ v (xf)(D°(xf)) = {0}, for all 
xj £ Pj, and therefore 

(A* - A*)| P/ = 0, (^ - ^)\p f = 0, (z/ r )| P/ = 0, 
for every i,j and r. This implies that (see (|70|)) 

(P(^))|p, = « " A^Zo))|p r (76) 

Finally, using (fT3|), (fr"j| ) and (|76"|), we conclude that 

[*(P(P^))^ ]| P/ = 0, [*(P(P^))r?°]| P/ = 1. 



Remark 14 If TPf n #^ v ((TP/)°) = {0} and jj : Pj — > J X *P is the canonical inclusion, then the 
couple (jj(o) v ), ff{fj)) defines a cosymplectic structure on Pf (see Remark |l2|) , and (F(E^ 7 ))\p f is 
just the evolution vector field of (h^ )|py with respect to the structure (jj(u) v ), ff{fj)). 



Remark 15 Let y' E be the vector field on E associated with another connection V' on tt : E —* P, 
and /ig 7 : Po — > M be the corresponding Hamiltonian function. Using Remark ||, we deduce that 

h?' = h v -iV 

is an extension to J l *E of the Hamiltonian function Kq' , where iV : J 1 *^ — > R is the real 
function on J 1 *P induced by the 7r-vertical vector field V = y' E — 3^ p . Thus, from Theorem [T^ 
and Proposition [irj (see Appendix we obtain that the evolution vector fields E^ and E^ v , 

of /i v and /i v ' with respect to the cosymplectic structures (u>^,fj) and (a) v ',r/) coincide. As a 
consequence, it follows that 

P(P^)=P(Pj',). 
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6.3 Dirac brackets and evolution of an observable 

Splitting TV J U E = D © #^ V (D°) allows us to introduce a Dirac bracket { , }^ , which can be 
defined as follows (see Appendix [B]). If F and G are C°°-differentiable real functions on U, then 

{F,G}% = ^(F(Xj),F(X'l)), 

where F : TjjJ^E — > D is the projector considered in Section 

A direct computation, using (p7\), proves that 

{F, G}% = {F, Gp + C^{F, X p }^{X a , G} v . (77) 

Moreover, we have 

Theorem 15 (i) The Dirac bracket { , }^ is a Poisson bracket and the Hamiltonian vector field 
of a C°° -differentiate real function F on U (with respect to { , }p ) is F(Xp). 

(ii) If g is an observable (that is, g is a C°° -differentiable real function on Pf), its evolution is 
given by the formula 

g = F(TZ^)\ Pf (G) + ({G,h^}S)\p f , 

where h v : U -> R and G : U -> R are arbitrary extensions to U of the Hamiltonian function 
/iq 7 : Po ~~ * K and of g : Pf — > R, respectively. 

(Hi) If F is a C°° -differentiable real function on U such that Xp(x) G j£^(D°)(x), for every 
x G U, then F is a Casimir function for the Dirac bracket, i.e., 

{F, G}p = 0, for every G G C°°(U). 

(iv) If F is a C°° -differentiable real function on U such that the restriction to Pf of Xp is tangent 
to P f , then ({F, G}J)\ Pf = ({F, G} v )\ Pf , for every G G C°°(U). 

( Proof) (i) Using Proposition [ll] (see Appendix |B|) , and the fact that D is a completely integrable 
distribution, we deduce that { , }^ is a Poisson bracket. Furthermore, from (|67| ) and (|77|), we 
obtain that 

F(Xj)(G) = {G, F}%, for F, G G C°°(U). 
(ii) If E^y = Vy ' -\-X^y is the evolution vector field of /i v with respect to the cosymplectic structure 



(u) v ,?7), then from Theorem 14 it follows that 

g = R(Ej)\p f (g). 

Thus, using the first part of this Theorem, we have that 

g = R(n V )\ Pf (G) + F(X^ v )\ Pf (G) = F(K*)\ Pf (G) + ({G, h v }^)\ Pf . 

(Hi) If F G C°°(U) and Xj(x) G # Av (D°)(x), for every x G U, then it is clear that F(Xj) = 0. 
Therefore, 

{F, G}% = -{G, F}~% = -F(Xj)(G) = 0, for every G G C°°(U). 

(iv) The condition (Xp)\ Pf G %(Pf) implies that (F(Xp))\ Pf = (Xp)\ Pf , and consequently 

({F,G}J)\ Pf = -( X J)\ Pf (G) = ({F,G}*)\ Pf . 
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From Theorem 15, we deduce that 

{X a , F}p = 0, 

for every a and F G C°°(U). In other words, the second class constraints {X a } are Casimir 
functions for the Dirac bracket { , }^ . Using Theorem 15, we also obtain that 



({fiV}?)|p, = ({CSVDIp,, {&\F}J)\ Pf = &,F}V)\ Pf , 

for every t and u, 1 < t < kof and 1 < u < kf — kof. 

Remark 16 If TPf n #7v(T(P/) ) = {0} and : P/ — > J^*E is the canonical inclusion, we can 
consider the cosymplectic structure (jj(u)) v ), J*f(fj)) on Pf whose Reeb vector field is P(^, v )|p / (see 
Remark ||). Furthermore, if X G £(£/) and G G C°°(f7), 

(i(P(X G v ))^ v )(P(X)) = (dG-^ v (G)r/)(P(X)) -cD v (Q(X G v ),P(X)). 

Thus, from (|67|), we have that 

(i(P(X G v )p v )(P(X)) = dG(F(X)) -P(7l v )(G)r?(P(X)) 

and, since P(X G )|p. is tangent to Pj, it follows that P(X G 7 )|p / is the Hamiltonian vector field of 
g = G\p f with respect to the cosymplectic structure (jj(a)) v ), J*t(fj))- Therefore, if on Pf (resp. 
J U E) we consider the Poisson structure induced by the cosymplectic structure (jj(cD) v ), J*f(fj)) 
(resp. by the Dirac bracket { , }p) then the canonical inclusion jf is a Poisson morphism. 



Remark 17 If V' is another connection on ir : E — > B, then we can obtain the same splitting of 
Tu(J u E) and the same projectors P : T v J l *E -» D and Q : Tc/J 1 *^ -» # Rv ,(D°) = # Av (D°). 
As a consequence, using (77) and Theorem Il2|, we conclude that the Dirac brackets coincide, i.e., 



{ , }d V 



{ , }J- 



7 Examples 



In this Section, we illustrate some aspects of the theory with two examples. In the first, we treat 
the particular case where the base manifold B is IR and the submersion tt is a trivial fibration of E 
on R and, as a consequence, we recover some results obtained in [ffH, |3l] . In the second example, 
we consider the special case of time-dependent Lagrangians which are affine on the velocities. 



7.1 Example 1 



Suppose that B = M and that 7r : E — > M is a trivial fibration; that is, Q is a differentiable manifold 
of dimension n, E = M x Q and tt : £7 = ix Q — > R is the canonical projection over the first factor. 
Then, the evolution phase space J l E of the system can be identified with the product manifold 
M. x TQ in such a way that the fibration tt 1 : J X E — > E is the map tq = Id x tq : E x TQ ^IxQ, 
where tq : TQ — > Q is the canonical projection (see f42|| ). 
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Now, let t be the usual coordinate on E and w the volume form given by w = dt. If £ is a 
Lagrangian density on J l E ^ MxTQ, we have that C = £n, with rj = (vf 1 )*^) and £ : RxTQ -> E 
the Lagrangian function. 

The Lagrangian function and the Liouville vector field C of TQ allow us to introduce the 
Lagrangian energy Eg : E x TQ — > M as the real function on M. x TQ defined by 

E c = C(£) - £. 

If (t, q u ,v u ) are fibered coordinates onMx TQ, we have that 

T-1 v d£ 

Ec = V d^~ £ - 

Furthermore, the vector field 3^ = §i (resp. y = Jj) on E x Q (resp. E x TQ) induces a 
distinguished connection V (resp. V) on the fibration tt : E x Q — > E (resp. tq : E x TQ — > 
E).) Thus, if Mj is the ith-Lagrangian dynamical constraint submanifold, one can give several 
descriptions of Mj as the zero set of the ith-generation dynamical Lagrangian constraints {C } 
(see Theorem |j). In particular, using Theorem |], we deduce that 

C W = i(zM)( V - j c ), for G X ± (M J ), 

where 7£ is the 1-form on E x TQ given by 7£ = £(Jt)Q,c, and fi£ is the Poincare-Cartan 2- form. 
This description of the constraints £W was obtained in [10]. On its turn, the description of the 



ith-generation of SODE Lagrangian constraints is given in Theorem 10 



Next, we will assume that the Lagrangian system (J l E,VLc) is almost regular. 

Note that in the particular case where the fibration tt is trivial, the restricted momentum dual 
bundle associated with tt can be identified with the product manifold E x T*Q. In addition, the 
map 

J U E x R = (M x T*Q) x E — ► T*E = T*(R x Q) 
((A, a),/i) — ► a + fj,dt(X) 

is a diffeomorphism. Moreover, in [|n]] it was proved that the Lagrangian energy Ejr is TCq- 
projectable onto a real function Kq : V — > E. Using this function, a diffeomorphism ho '■ V Q 
E x T*Q -^PC T*£ = (E x T*Q) x E can be defined as follows 

foo(ai) = (5, — fyf(x)), for every x £ V. 

A direct computation shows that ho = ho = fin , where fiQ : V — > T 7 is the restriction to "P of the 
canonical projection /i : T*(M x Q) — > E x T*Q (see Section |4.1| ). Therefore, the Hamilton-Cartan 
forms ®h and O/j are given by = (Jo ° ho)*@ and = (Jo ° fy))*^> © and being the 
Liouville 1-form and the Liouville 2-form, respectively, of T*(M x Q). 

Now, if y is a vector field on V which is TQ-projectable onto the vector field 3^ = Jj, then y 
induces a connection V in the fibration Tq : V — » E. Using this connection, one can give several 
descriptions of the ith-Hamiltonian dynamical constraint submanifold Pj as the zero set of the ith- 
generation dynamical Hamiltonian constraints (see Theorem 0). In particular, from Theorem 
0, it follows that 

£0 = - 7ho ), for Z« G X ± (P,), 

where 77 = (fQ)*dt and ^h is the 1-form on V given by 7/^ = i(y)Qh - This description of the 
constraints £W was obtained in fpTof] - 
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Moreover, it is clear that /ig 7 is just the Hamiltonian function Kq on V associated with the 
Lagrangian system, the connection V on 7r : E — * M and the 1-form w = dt. In addition, if 
is an extension to M. x T*Q of the Hamiltonian function , then we deduce 
that V and u) v are the Liouville 1-form and the Liouville 2-form of T*Q, respectively, that 
K v = §- v and that 

9 h v = G v + h^dt, n h v = w v + dh v A dt. 

As a consequence, the bracket { , } v onIx T*Q is the usual Poisson bracket on K x T*Q (see 
Section 16.11). 



Finally, suppose that {A^,} is a maximal set of local independent second class constraints for 
the final Hamiltonian dynamical constraint submanifold Pf, and denote by C a p the real function 
defined by 

C a p = C a/3 + — — — = {X a , Xp) + — — for every a and (5. 

Then, since the matrix (C a p) is regular, we deduce that the matrix (C a p) also is regular. In fact, if 
(C a P) is the inverse matrix of (C a p) and 

r af3 _ na/3 , part nfiv 

L ~ L +L L dt dt 

then (C a @) is the inverse matrix of (C Q p) (note that the skew-symmetric character of the matrix 

- - / / dX i QX i 

(C a P) implies that C 7 v — = 0). Thus, if { , }^ is the corresponding Dirac bracket on 



K x T*Q, using fl77Q , we conclude that 

{F,G}% = {F,G} V -C aP {Xp,G}^{X a ,F} v 

-CapC a 'l3'{Xl3,G} V '{Xpt,F} V ~g^ Qt ' 

for F,G € C°°(M x T*Q). This expression for the Dirac bracket { , }^ was obtained in [31] 



7.2 Example 2 

Let 7 be a 1-form on the configuration space E and £ : J 1 E — > R be the Lagrangian function 
defined by 

£(jt <l>) = i(jt <t>) = im)), for j/ feJ'E. 

If we take fibered coordinates (t, q l , Vj) on J l E such that 7 = 7j(i, g)<i(f + 70 (i, then 

7(*, = 7i(MK+7o (*,?)• 

Thus £ = 7 is an affine Lagrangian on the velocities. 

A direct computation shows that ©£ = (7r 1 )*7 and hence = —d&£ = (vr 1 )*(— dj). We also 
obtain that J-C = 7 o 7T 1 . Therefore, V is an embedded submanifold of T*E which is diffeomorphic 
to E. The mapping J-Cq : J^E — > P = 7(-E) may be viewed as the composition 

J X E ^ E 7(E) = V. 

Since 7 : E —* 7(-E') is a diffeomorphism and 7T 1 is a surjective submersion with connected fibers, 
we deduce that J 7 Co '■ J 1 ^ — » "P is also a surjective submersion with connected fibers. Moreover, 
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—l 



it is easy to prove that TCo (FCq{x)) = J 7 C 1 (J r Co(x)) = (ir 1 ) 1 (ir 1 (x)), for every x € J l E. As 
a consequence, the Lagrangian is almost-regular. Note that the diagram 



V = FC {J l E) 




^E 



is commutative, that the mappings 7 : E — > V and /iq : V — * V are diffeomorphisms and that 
(Mo 7)*(^h ) = ~~ ^7> an d (a*o 7 )*( r ?°) = 7r*(ro). Thus, the submanifold V can be identified with 
E and, with this identification, the mapping J- Co : J^E — > "P is the submersion -zr 1 : J 1 ^ — > E 1 and 
the 2-form f2/j (resp. the 1-form rp) on P is the 2-form — dj (resp. the 1— form 7r*(dt)). Taking 
these identifications into account, the constrained Hamilton equations on E are 



i(r)(d 7 ) = o, 



i(Y)(ir*(dt)) 



1. 



(78) 



Next, we will assume that the couple (c?7, 7r*(df)) is a cosymplectic structure on E with Reeb vector 



field ft. 
and that 



. Then, a direct computation shows that the matrix (7^,- = -—4- — — — )i<i 7 < n is regular 



7 5t + 2 7 dt J dqj' 



(79) 



where (7* 



)l<i,j<n 



is the inverse matrix of (7 



ij )l<i,j<n- 



Furthermore, it is clear that the final constraint submanifold for the Hamiltonian problem is 



Po = V and that there exists a unique solution Yh of the constrained Hamiltonian equations (78); 
namely, the Reeb vector field TZj. 

Moreover, if Xc is a vector field on J l E which projects via J- Co onto 7£ 7 = Y^ , then Xc is a 
solution of the dynamical Lagrangian equations on J l E. In addition, the corresponding submanifold 
S of J l E (see Section |4.5|) is IZj(E) and the unique vector field X^ = Tc on S satisfying 



i(Tc)^c\s = 0, i(T c )v\s = l, 



J(r c )\s = o 



is T c = (ft 7 

results of Section 4.5, we deduce that S can be locally described as follows 



) c \s, where (7£ 7 ) c denotes the complete lift of 7£ 7 to TE (see |32||). From (|79|) and the 



S={(t,q j ,v j ) G J 1 ^ 



77"' 



dlo 
dq l 



dt 



Vi}. 



Another way of finding Euler-Lagrange vector fields is by applying the procedure described in 
Section ||. As there is solution of the dynamical Lagrangian equations ([!]) everywhere in J l E, 
we look for the submanifold of J l E where solutions satisfying the SODE condition exist. The 
corresponding compatibility conditions (lst-generation SODE Lagrangian constraints) are given 
by Prop. 0. In this case, these constraints are 



0«i 



yj ryW 



1 dq i 



ot 



Hence S\ 



S. The corresponding Euler-Lagrange vector fields are 

'#70 _ <9tA d 
dq l dt 



d 1 

1 — 7 

dt 2 1 



Kl 



dqi dvi 
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The tangency conditions ^£((f> )\si = determine all the coefficients , then we obtain the 
unique solution Tc\s-i = C^"y) c |si- Thus, Sf = Si = S (see the comments made in Section 



Now, let V be a connection on tt : E —* B, and y E be the associated vector field on E such 
that TT*(w)(y E ) = 1. 



Using ( |5l| ) and (|5^), we obtain that 

h% = -i(y E ), 

fly : Pq = E — > R being the Hamiltonian function. 

Moreover, a set of local independent constraint functions defining Pq = V as a submanifold of 
J X *E is given by 

d° } = Pi -7. (*.?), l<i<n. (80) 



Prom (||) and (|(| it follows that 

I Ipo(-e)-^ -^J -7ii (81) 

and, hence we conclude that the matrix } V |p )i<« i i<« ^ s regular. Thus, all the constraints 

{£i }i<*<n are second class and TPo H #^ v ((TPo)°) = {0}. Furthermore, with the identification 
P = V = E, the cosymplectic structure on E 1 induced by the couple (u3 v ,r/) (see Remark 12) is 

(_d 7 + d( 7 (y B )) A ir*(dt),Tr*(dt)). 

Finally, if { , j/j 7 is the Dirac bracket on J U E, then using (|6^) , (|77| ) and fl8l]), we deduce that 

v _ fd£dG _ 8F_dG\ _ ^ ( dF t d 7i dF\ ( dG , d 7i 



1 ' i; ^dq l dpi dpidq 1 ) ^ \dq % dq k dpk) \dqi dq l dp, 



for F,Ge C^iJ^E). 
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Appendices 

A Precosymplectic vector spaces 

Let V be a real vector space of finite dimension, and suppose that u> is a 2-form on V and r/ £ V*. 
We can consider the linear map \>(n,u) '■ V ^ V* given by 

b(r),Lu){ v ) = i(v)u> + (i(v)r])r), for every v G V. 

If S is a subspace of V, we will denote by S° the annihilator of S; that is, 

S° = { a e V*/i(u)a = 0, VnG S}, 
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and, in the same way, for a subspace W of V*, 

W° = {u£ V/i(u)a = 0, VqG W}. 
Then, if K is a subspace of V we define the orthogonal subspace of K with respect to r/ and cu as 
^ = (W^))° = i v e t7(i(«)<" " (i^MI* = 0}. (82) 

A precosymplectic structure on a real vector space on V is a triad (rj,u,TZ), where r] £ V* , uj 
is a 2-form on V, H G V and 

i{K)uj = 0, i(7e)r? = 1. 
Note that, in this case, the rank of bfau) is odd. Moreover, we have 

Lemma 12 Let (rj,u,TZ) be a precosymplectic structure on a real vector space. Then: 

(i) V = ker ui n ker r\. 

(ii) If K is a subspace of V , then dim K 1 - = dim V — dim K + dim(y _L n K) . 
(Hi) IfK and K' are subspaces ofV, K' C K then K L C (K') L . 

( Proof) (i) It is clear that kerw fl kern C V- 1 -. Conversely, if u G F 1 " then 

= (i(u)w - (i(v)n)ri)(K) = -i(v)n, 
which implies that i(v)u = 0, that is, v G kerw n ker 77. 
(ii) Consider the linear map 

^,u,)|K = # -» V"*, u£K -» b(„^)(«) = + (i(u)n)ri G V*. 

A direct computation shows that ker(b( r? w )|^-) = fl .fT. Moreover, we have 

= W*) = ((b M (K))°)° = (K^)°. 

Therefore, 

dim if - dim(l/ ± n if) = dimF - dimK ± . 

(Hi) It follows from ( jB^) . ■ 

Now, let a 2-form on a real vector space V of finite dimension and rj G V* . If 7£ is a vector 
of V such that r/(7£) = 1, we can consider the 2-form ujh defined by 

un = ^ - V A 77?. (83) 

where 77^. G is given by 

77 e = i(K)(n). (84) 
It is clear that the triad (n,uJii,7Z) is a precosymplectic structure on V. 

If X is a subspace of V, we will denote by ^ the orthogonal subspace of X with respect 
to uj-ji and 77. 
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Lemma 13 Suppose that TZ and TZ' are vectors of V such that rj(TZ) = n(TZ') = 1. Then: 

(i) V,\ , = v, 1 - 



Kfl) i w n', n )' 
(ii) If K is a subspace of V, we have 



dim K tl , „\ = dim K ( , , „\ . 



( Proof ) (i) We must prove that ker Wft n ker rj = ker u-r; fl ker rj (see the first part of Lemma 12) . 

Prom (|8^) and (f3), we deduce that 

u-R> = + 7] A (7^ - 7^/). 
Thus, if v € kerwft n ker rj then i(v)ui-jn = — i(v)(^-ji — 7ft' )/?. In particular, 

= i(v)uK> = -i(u)(7ft -7w) 
which implies that i(v)u>n' = 0. 

This shows that kerwft D ker?] C kercjft' n kerr?. The converse is proved in a similar way. 



(ii) It follows from (i) and the second part of Lemma 12 



Next, we will exhibit the definition and some properties of cosymplectic structures on a real 
vector space (for more information, see [§, |], |3l|). 

Let V be a real vector space of dimension 2n + 1. 

A cosymplectic structure on 1/ is a couple (rj,u), where n E V*, u is a 2- form on and the 
linear map \>(n,uA '■ V ^ V* is a. linear isomorphism or, equivalently, rj Aw" = r\ Auj A . . S n . . . Auj ^ 0. 
If (V, n, uj) is a cosymplectic vector space then TZ = ^J^^ (rj) is called the Reeb vector of V. TZ is 
the unique vector of V which satisfies the conditions 

i(TZ)n = 1, i(TZ)u = 0. 

In particular the triad (rj, ui, TZ) is a precosymplectic structure on V . Furthermore, using Lemma 
|l2^ , it follows that 



Lemma 14 t 31] Let (rj,ui) be a cosymplectic structure on a real vector space V . Then: 



1. v 1 - = {0}. 

2. If K is a subspace of V , we have that dim K = dim V — dim K. 

3. If K is a subspace of V and K n K ± = {0} then V = K®K L . 

We also get 

Lemma 15 Let (rj, ui) be a cosymplectic structure on a real vector space V . Suppose that K is a 
subspace of V such that K n K = {0} and that there exists TZk G K satisfying 

i(7Z K )u K = 0, i(TZ K )rjK = 1, (85) 

where ujk = i^\kxK °> n d r/x = v\k- Then the couple (tjk-i^k) is a cosymplectic structure on K and 
TZk is the Reeb vector of K. 
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( Proof) Let ^Oq K ,u K ) '■ K — > K* be the linear map induced by the 2-form ujk and the 1-form rjx- 
If X G ker br VKjUlK ) then we deduce that i{X)ujk = 0, i(X)rjx = 0. Now, since K n K 1 - = {0}, it is 
easy to prove that 

ker uj K n ker n K = {0} (86) 
and thus X = 0. Therefore, the map ^(r/ K ,uj K ) is a linear isomorphism. 

Furthermore, if X G keru>K then, from (p5j), it follows X — t]k{X)TZk G kerwx n ker t)k and 
consequently, from (|86|), we obtain that X = T]k(X)1Zjc. Thus, kerux = (TZk), which implies that 
the dimension of K is odd. ■ 



B Poisson and cosymplectic structures and Dirac brackets 

A Poisson structure on a differentiable manifold M is a bracket of functions { , } : C°°{M) x 
C°°(M) — > C°°(M) which satisfies the following properties: 

• Skew-symmetry: {F, G} = — {G, F}. 

• Leibniz rule: {FF', G} = F{F', G} + F'{F, G}. 

• Jacobi identity: {F, {G, H}} + {G, {H, F}} + {H{F, G}} = 0. 

If { , } is a Poisson bracket and F G C°°(M) one can define the Hamiltonian vector field Xp 
associated with F as follows 

X F (G) = {G, F}, for every G G C°°(M). 

In addition, one also can define a 2- vector A characterized by the condition 

A(dF, dG) = {F, G}, for F,G e C°°(M), (87) 

and it follows that [A, A] = 0, where [ , ] is the Schouten-Nijenhuis bracket. 

Conversely, if A is a 2-vector on M such that [A, A] = and { , } is the bracket of functions 



given by (87), then { , } is a Poisson bracket (for more details, see j33l §j§). 



Next, we show that a cosymplectic structure on manifold M induces a Poisson bracket on 
C°°(M) (see [§, [|]). We need some results about cosymplectic vector spaces. 

Let (V, 7j, u) be a cosymplectic vector space with Reeb vector TZ. If a G V*, we can consider the 
vectors X a and E a defined by 

These vectors are characterized by the following conditions 

i(X a )u = a — a(TZ)r), i(X a )rj = 0, 
i(E a )u> = a — a(TZ)r), i(E a )n = 1. 

Note that E a =K + X a . 

The cosymplectic structure (rj, u>) also allows us to introduce a 2-vector A : V* x V* — > M on V 
given by 

k{a,f3)=u{X a ,Xp)=u{E a ,E p ), fora,/3GV*. 
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The 2-vector A induces a linear map #a ■ V* — ► V denned by /3(#a(o0) = A( a >/^)> f° r &,f3 £ V*. 
We have that 

t , ( f? ,o;)(#A(a)) = a - a{K)r). 
In particular, 4f\(V*) = ker 77, ker #a = ( r ?) and the map #a : (7^)° ~~ > ker 77 is a linear isomorphism. 

If if is a subspace of V which satisfies certain conditions, then the map #a can be used in order 
to obtain the orthogonal subspace of if with respect to (77, oS). 

Lemma 16 If if is a subspace ofV and there exists TZk £ if such that 

(i(R K )u)\ K = 0, i{n K )n = l (88) 

then if- 1 = #a(K°). 

( Proof) Suppose that a £ if . Since rj(jf\{a)) = 0' ^ follows that 

(i(# A («))w-»7(#A(a))»/)|jr = «#A(a)H|x = (b ( ^)(# A (a))k 

= (a-a(ft)»j)|tf = (-a(ft)fj)|jf. 

Therefore, from (pq), we deduce that 

= (-a(7ty7)(ftjc) = -a{K) 
and (i(# A (a))w - n(# K (a))rj)\ K = 0, i.e, # A (a) € if- 1 . We have proved that #a(K°) C if- 1 . 
Now, using Lemma |l4] and the fact that if n (77) = {0}, we obtain that 

dim# A (#°) = dim if = dim if- 1 . 



Moreover, a direct computation shows the following. 

Proposition 9 Let (n,u>) be a cosymplectic structure on a real vector space V and 1Z be the Reeb 
vector ofV. Suppose that a is a 1-form on V and that u> is the 2-form on V defined by Q = tu+aAn. 
Then: 

1. The couple (r),ui) is a cosymplectic structure on V with the Reeb vector 1Z given by 7Z = E a , 
where E a = + (1 - a(K))r}). 

2. If (3 is a 1-form on V and Xp (resp. Xp ) is the vector ofV defined by Xp = \)7^A(3 — fiiJVjrj) 
(resp. Xp = } (/9 - P(K) V )) then Xp = Xp. 

3. If A (resp. A) is the 2-vector on V induced by the structure (rj, uS) (resp. (r/,cD)) then A = A. 

Now, let M be an odd-dimensional differentiable manifold, and (r/, u) a cosymplectic structure on 
M; that is, 77 is a closed 1-form, to is a closed 2-form and the couple (r] x ,uj x ) is a cosymplectic 
structure on the vector space T X M, for every x £ M. Then there exists a unique vector field 
1Z on M, the Reeb vector Eeld, satisfying the conditions i(JZ)uj = and i(TZ)n = 1. Moreover, 
using the above constructions on cosymplectic vector spaces, one can take the 2-vector A on M 
induced by the couple (r},oj), and it follows that A defines a Poisson structure on M (sec ]|, [|]). In 
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addition, if F € C°°(M) we also can consider the vector fields Xp and Ep on M characterized by 
the conditions 

i{X F )uj = dF-K(F)r], i{X F )n = 0, 
i{E F )uj = dF-K(F)ri, i{E F )rj = 1. 

Xf is just the Hamiltonian vector field of F with respect to the Poisson structure A. Ef is called 
the evolution vector field of F with respect to the cosymplectic structure (t),u). We have that 
E F = 11 + X F . 

Furthermore, from Proposition |9], we deduce the following result. 

Proposition 10 Let (rj, u) be a cosymplectic structure on a differentiable manifold M and 1Z 
be the Reeb vector field of M. Suppose that F is C 00 -differentiable real function on M and that uj 
is the 2-form on M defined by uj = uj + dF A r/. Then: 

1. The couple (rj,Q) is a cosymplectic structure on M with Reeb vector field 1Z given by 

a = E F , 

where Ep is the evolution vector field of F with respect to the cosymplectic structure (ji,u)). 

2. If G is a C°° -differentiable real function on M and Xq (resp. Xq) is the Hamiltonian vector 
field of G with respect to the structure (t],u) (resp. (rj,ui)) then Xq = Xq. 

3. If A (resp. A) is the Poisson 2-vector on M associated with the cosymplectic structure (r],Lu) 
(resp. (ri,u>)) then A = A. 

Next, assume that (r/, uj) is a cosymplectic structure on a manifold M and that D is a distribution 
on M of dimension d satisfying 

TM = J D0# A ( J D°), (89) 

A being the Poisson 2-vector associated with the structure {ji,uj). The splitting ( |89[ ) allows us to 
introduce two projectors 

P : TM -» D, Q : TM -> # A (D°) 
in such a way that X = FX + QX, for every X G TM. 

Moreover, one can define a Dirac bracket of functions { , }d as follows 

{F,G} d =u;(¥(Xp),¥(Xg)), 

for F,G£ C OQ (M), where Xp and Xq are the Hamiltonian vector fields of F and G with respect 
to the cosymplectic structure (77, uj). It is clear that the bracket { , }p> is skew-symmetric and that 
it satisfies the Leibniz rule. 

Moreover, since (?/(#a(«)) = for every a, we deduce that rj(P(TZ)) = 1. Thus, it is easy to 
prove that D{x) = D(x) © (P(TZ)(x)) , for every x € M, where D{x) = D(x) Hker n x . Therefore, the 
assignment 

x E M — > D(x) C T X M 
defines a distribution D on M of dimension d — 1. In addition, we have 

Proposition 11 J^J The distribution D is completely integrable if and only if { , }o is a Poisson 
bracket. 

Remark 18 If the distribution D is completely integrable then it is clear that the distribution D 
is also completely integrable. 
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C Auxiliary results 

Taking into account the characterization of the submanifolds Cj given in Proposition |], we can 
prove the following results (remember that Cq = F): 

Lemma 17 For every x G Cq, ker Q, x n ker r\ x C ker u x Pi ker rj x = T x Cq. 

( Proof ) The last equality follows from Lemma |l^. 
For the first inclusion, if X G ker Q x n ker 77 x then 
i(X)w x = »(X)J2 X - i(X)(n x A 7x ) = i{X) lx f\rj x = i(X) i(y)tl x A r/. 
and thus X G T x Cq. 

Lemma 18 For every x G C\, ker H ker 77^ = ker cj x n ker = T x Cq. 

( Proof ) If X G T^C = ker u x n ker % then 

= 7(X)cj x + i(A")(7/ a A 7x ) = -[t(X)7 x ]7fe 

but as 1 G Ci, there exists Y G T x Co such that i(Y)Q x = and i(Y)r/ x = 1; therefore the 
contraction of the above equality with this Y leads to i(X)j x = 0, and hence i(X)Q x = 0. ■ 

Lemma 19 For every Z,Y G ker cj n ker 77, we get 

1. i(Y)i(Z)d 7 | Cl =0. 

2. [Z,Y]| Cl G (ker cj n ker r?)| Cl . 

( Proof) 

1. As Z, Y G ker cj n ker 77, we obtain that i(Y)i(Z)fi = i(Y)i(Z)cj + i(Y)i{Z)(n A 7 ) = 0. 
Therefore 

i(Y)i(z)d 7 = i(y)»(z)dt(y)n = i(y)t(^)[L(y)n] 

= L(y) - i(y) z])n - y]) *(z)o 

= i([^,Z])i(Y)0-i([^Y])i(Z)0 

However, as a consequence of the above Lemma, since Z\c 1 ,Y\c 1 G (ker f2 n ker 77) [cr a , hence 
i(Y)i(Z)d 7 | Cl =0. 

2. Observe that as Z, Y G ker cj D ker 77, and 77 is a closed form, 

i([Z, Y])r, = L(Z) i(Y) V - i(Y) L(Z) ?? = - i(Y) i(Z)d V - i(Y)di{Z) V = , 
then [Z, Y] G ker 77. Moreover, 

i([Z, Y])cj = L(Z) i(Y)cj - i(Y) L(Z)cj = - i(Y) i(Z)dcj - i(Y)d?;(Z)cj = - i(Y) i(Z)du. 



= -i{y)i{X)fl x Ar lx = 
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Now, for every X € X(Cq), as du: = df2 — d(n A 7) = rj A dj, we have 

i(X) i(Y) i{Z)du = i(X) i(Y) i(Z)(rj A d 7 ) = [i(X)rj] i(Y) i{Z)d 1 . 
Therefore, for every x G C\, taking into account the first item we obtain 

i(X x )i(Y x )i(Z x )(dj) x = [i(X x ) Vx ]i(Y x )i(Z x )(d 7 ) x = 
and the result follows. 
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